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Abstract. This paper is devoted to the study of mean-field limit for systems of in- 
distinguables particles undergoing collision processes. As formulated by Kac [22] this 
limit is based on the chaos propagation, and we (1) prove and quantify this property 
for Boltzmann collision processes with unbounded collision rates (hard spheres or long- 
range interactions), (2) prove and quantify this property uniformly in time. This yields 
the first chaos propagation result for the spatially homogeneous Boltzmann equation for 
true (without cut-off) Maxwell molecules whose "Master equation" shares similarities 
with the one of a Levy process and the first quantitative chaos propagation result for 
the spatially homogeneous Boltzmann equation for hard spheres (improvement of the 
convergence result of Sznitman [40]). Moreover our chaos propagation results are the 
first uniform in time ones for Boltzmann collision processes (to our knowledge), which 
partly answers the important question raised by Kac of relating the long-time behavior 
of a particle system with the one of its mean-field limit, and we provide as a surprising 
application a new proof of the well-known result of gaussian limit of rescaled marginals 
of uniform measure on the JV-dimensional sphere as N goes to infinity (more applications 
will be provided in a forthcoming work). Our results are based on a new method which 
reduces the question of chaos propagation to the one of proving a purely functional es- 
timate on some generator operators {consistency estimate) together with fine stability 
estimates on the flow of the limiting non- linear equation {stability estimates). 
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1. Introduction and main results 

1.1. The Boltzmann equation. The Boltzmann equation (Cf. [10] and [11]) describes 
the behavior of a dilute gas when the only interactions taken into account are binary 
collisions. It writes 

(1.1) %+v-V x f = Q(fJ) 

where Q = Q(f, f) is the bilinear Boltzmann collision operator acting only on the velocity 
variable. 

In the case when the distribution function is assumed to be independent on the position 
x, we obtain the so-called spatially homogeneous Boltzmann equation, which reads 

(1.2) ^(t,v) = Q(fJ)(t,v), veR d , t>0, 

where d > 2 is the dimension. 

Let us now focus on the collision operator Q. It is defined by the bilinear symmetrized 
form 

(1.3) Q{gJ){v) = \ I B(\v-v*\,cos6)(g'J' + g'fl-g*f-gf*)dv*do-, 

where we have used the shorthands / = f(v), f = f(v'), g* = <?(«*) and «j£ = g{v'j). 
Moreover, v' and v[ are parametrized by 

(1-4) v = — ^— + 2 ' a, < = — 2 ' a 

Finally, 6 G [0, ir] is the deviation angle between v' — v'^ and v — defined by cos 8 = a -u, 
u = v — u*, u = u/\u\, and B is the Boltzmann collision kernel determined by physics 
(related to the cross-section S(u — v#, a) by the formula B = \v — v* \ S). 

Boltzmann's collision operator has the fundamental properties of conserving mass, mo- 
mentum and energy 



(1.5) / Q(f,f)<f>(v)dv = 0, <t>(v) = l,v,\v 
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and satisfying the so-called Boltzmann's H theorem 

We shall consider collision kernels B = T(\v — u*|) 6(cos 9) (with T, b nonnegative func- 
tions). Typical physical interesting kernels are in dimension 3: 

• (HS): the hard spheres collision kernel B{\v — cos 6) = est \v — 

• collision kernels deriving from interaction potentials V(r) = cstr _s , s > 2: T(z) = 
|z| 7 with 7 = (s — 4)/s, and b is L 1 apart from 9 ~ 0, where 

b(cos9) ~6»~o Cf, # _ ^~ 1 '~ z/ with v = 2/s (see [10]) including in particular: 

(tMM): the true Maxwell molecules collision kernel when 7 = and v = 1/2; 

• (GMM): Grad's cutoff Maxwell molecules when B = 1. 
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1.2. Deriving the Boltzmann equation from many-particle systems. The question 
of deriving the Boltzmann equation from particles systems (interacting via Newton's laws) 
is a famous problem, related to the so-called 6-th Hilbert problem mentionned by Hilbert 
at the International Congress of Mathematics at Paris in 1900. 

At least at the formal level, the correct limiting procedure has been identified by 
Grad [18] in the late fourties (see also [9] for mathematical formulation of the open ques- 
tion): it is now called the Boltzmann- Grad or low density limit. However the original 
question of Hilbert remains largely open, in spite of a striking breakthrough due to Lan- 
ford [25], who proved the limit for short times. The tremendous difficulty underlying 
this limit is the irreversibility of the Boltzmann equation, whereas the particle system 
interacting via Newton's laws is a reversible Hamiltonian system. 

In 1954-1955, Kac [22] proposed a simplified problem in order to make mathematical 
progress on the question: start from the Markov process corresponding to collisions only, 
and try to prove the limit towards the spatially homogeneous Boltzmann equation. Going 
back to the idea of Boltzmann of "stosszahlansatz" (molecular chaos), he formulated the 
by now standard notion of chaos propagation. 

Let us first define the key notion of chaoticity for a sequence (f N )N>i of probabilities 
on E , where E is some given Polish space (and we will take E = M. d in the applications): 
roughly speaking it means that 

f N ~ f® N when N^oo 

for some given one-particle probability / on E. It was clear since Boltzmann that in the 
case when the joint probability density / of the TV-particles system is tensorized into N 
copies f® N of a 1-particle probability density, the latter would satisfy the limiting Boltz- 
mann equation. Then Kac made the key remark that although in general coupling between 
a finite number of particles prevents any possibility of propagation of the "tensorization" 
property, the weaker property of chaoticity can be propagated (hopefully!) in the correct 
scaling limit. The application example of [22] was a simplified one-dimensional collision 
model inspired from the spatially homogeneous Boltzmann equation. 

The framework set by Kac is our starting point. Let us emphasize that the limit 
performed in this setting is different from the Boltzmann-Grad limit. It is in fact a mean- 
field limit. This limiting procedure is most well-known for deriving Vlasov-like equations. 
In a companion paper |34] we shall develop systematically our new functional approach 
for Vlasov equations, McKean-Vlasov equations, and granular gases Boltzmann collision 
models. 

1.3. Goals, existing results and method. Our goal in this paper is to prove (and 
set up a general robust method for proving) chaos propagation with quantitative rate in 
terms of the number of particles N and of the final time of observation T. Let us explain 
briefly what it means. The original formulation of Kac [22] of chaoticity is: a sequence 
f N E P sym (E N ) of symmetric probabilities on E N is /-chaotic, for a given probability 
/ G P(E), if for any £eW and any tp E C b {E)® 1 there holds 

!*->-</•«.*) 

which amounts to the weak convergence of any marginals (see also [6] for another stronger 
notion of "entropic" chaoticity). Here we will deal with quantified chaoticity, in the sense 
that we measure precisely the rate of convergence in the above limit. Namely, we say that 
f N is /-chaotic with rate e(N), where e(N) — > when N — > oo (typically s(N) = N~ r , 
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r > 0), if for some normed space of smooth functions F C Cb(E) (to be precised) and for 
any t G N* there exists Kg G (0,oo) such that for any ip G J 7 ®^, \\<p\\f < 1, there holds 

(1-6) \(ll e [f N ]-f^^)\<K e e(N), 

where IT^ [f N ] stands for the ^-th marginal of / . 

Now, considering a sequence of densities of a ^-particles system f N G C([0, oo); P sym {E N )) 
and a 1-particle density of the expected mean field limit / G C([0, oo); P(E)), we say that 
there is propagation of chaos on some time interval [0, T] if the /o-chaoticity of the initial 
family /q implies the /t-chaoticity of the family f^ for any time t G [0, T]. 

Moreover one can roughly classify the different questions around chaos propagation into 
the following layers (in parenthesis, the corresponding probabilistic interpretation for the 
empirical measure, see below): 

(1) Proof of propagation in time of the convergence in the chaoticity definition (prop- 
agation of the law of large numbers along time). 

(2) Same result with a rate e(N) as above (estimates of the rates in the law of large 
numbers, estimates on the size of fluctuations around the deterministic limit along 
time). 

(3) Proof of propagation in time of the convergence to a "universal behavior" around 
the deterministic limit (central limit theorem). See for instance |31| [38] for related 
results. 

(4) Proof of propagation in time of bounds of exponential type on the "rare" events 
far from chaoticity (large deviation estimates). 

(5) Estimations (2)-(4) can be made uniformly on intervals [0, T) with T finite or 
T = +oo. 

For Boltzmann collision processes, Kac |22j-|23| has proved the point (1) in the case of 
his baby one-dimensional model. The key point in his analysis is a clever combinatorial 
use of a semi-explicit form of the solution (Wild sums). It was generalized by McKean 
|32| to the Boltzmann collision operator but only for "Maxwell molecules with cutoff", 
i.e., roughly when the collision kernel B above is constant. In this case the combinatorial 
argument of Kac can be extended. Kac raised in |22] the question of proving chaos 
propagation in the case of hard spheres and more generally unbounded collision kernels, 
although his method seemed impossible to extend (no semi-explicit combinatorial formula 
of the solution exists in this case). 

In the seventies, Griinbaum |20] then proposed in a very compact and abstract paper 
another method for dealing with hard spheres, based on the Trotter-Kato formula for 
semigroups and a clever functional framework (partially remindful of the tools used for 
mean-field limit for McKean- Vlasov equations). Unfortunately this paper was incomplete 
for two reasons: (1) It was based on two "unproved assumptions on the Boltzmann flow" 
(page 328): (a) existence and uniqueness for measure solutions and (b) a smoothness 
assumption. Assumption (a) was indeed recently proved in |17] using Wasserstein metrics 
techniques and in [14] adapting the classical DiBlasio trick [13] , but concerning assumption 
(b) , although it was inspired by cutoff maxwell molecules (for which it is true) , it fails for 
hard spheres (cf. the counterexample built by Lu and Wennberg in [29] ) and is somehow 
"too rough" in this case. (2) A key part in the proof in this paper is the expansion of the 
a Hf" function, which is an a clever idea of Griinbaum (and the starting point for our idea 
of developing an abstract differential calculus in order to control fluctuations) — however 
it is again too rough and is adapted for cutoff Maxwell molecules but not hard spheres. 
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A completely different approach was undertaken by Sznitman in the eighties [JT] (see 
also Tanaka |43|), Starting from the observation that Griinbaum's proof was incomplete, 
he gave a full proof of chaos propagation for hard spheres. His work was based on: (1) 
a new uniqueness result for measures for the hard spheres Boltzmann equation (based 
on a probabilistic reasoning on an enlarged space of "trajectories"); (2) an idea already 
present in Griinbaum's approach: reduce by a combinatorial argument on symmetric 
probabilities the question of chaos propagation to a law of large numbers on measures; (3) 
a new compactness result at the level of the empirical measures; (4) the identification of 
the limit by an "abstract test function" construction showing that the (infinite particle) 
system has trajectories included in the chaotic ones. Hence the method of Sznitman proves 
convergence but does not provide any rate for chaoticity. Let us also emphasize that 
Graham and Meleard in [19] have obtained a rate of convergence (of order 1/y/N) on any 
bounded finite interval of the A-particles system to the deterministic Boltzmann dynamic 
in the case of Maxwell molecules under Grad's cut-off hypothesis, and that Fournier and 
Meleard in [151 116] have obtained the convergence of the Monte-Carlo approximation (with 
numerical cutoff) of the Boltzmann equation for true Maxwell molecules with a rate of 
convergence (depending on the numerical cutoff and on the number N of particles). 

Our starting point was Griinbaum's paper [20]. Our original goal was to construct a 
general and robust method able to deal with mixture of jump and diffusion processes, as 
it occurs for granular gases (see for this point the companion paper [33J ) . It turns out that 
it lead us to develop a new theory, inspiring from more recent tools such as the course of 
Lions on "Mean-field games" at College de France, and the master courses of Meleard [33] 
and Villani [35] on mean-field limits. One of the byproduct of our paper is that we make 
fully rigorous the original intuition of Griinbaum in order to prove chaos propagation for 
the Boltzmann velocities jump process associated to hard spheres contact interactions. 

As Griinbaum [20] we shall use a duality argument. We introduce the semigroup 
associated to the flow of the A-particle system and its "dual" semigroup. We also 
introduce the (nonlinear) semigroup associated to the meanfield dynamic (the expo- 
nent "NL" recalling that the limit semigroup is nonlinear in the most physics interesting 
cases) as well as T t °° the associated (linear) "pushforward" semigroup (see below for the 
definition). Then we will prove the above kind of convergence on the linear semigroups 
T t N and T t °°. 

The first step consists in defining a common functional framework in which the N- 
particles dynamic and the limit dynamic make sense so that we can compare them. Hence 
we work at the level of the "full" limit space P{P{E)) (see below). Then we shall identify 
the regularity required in order to prove the "consistency estimate" between the generators 
G N and G°° of the dual semigroups T t N andr t °°, and then prove a corresponding "stability 
estimate" at the level of the limiting semigroup Sj^ L . The latter crucial step shall lead 
us to introduce an abstract differential calculus for functions acting on measures endowed 
with various metrics. 

In terms of existing open questions, this paper solves two related problems. First it 
proves quantitative rates for chaos propagation for hard spheres and for (non cutoff) 
Maxwell molecules. These two results can be seen as two advances into proving chaos 
propagation for collision processes with unbounded kernels (in the two physically relevant 
"orthogonal" directions: unboundedness is either due to growth at large velocities, or to 
grazing collisions). Second we provide the first uniform in time chaos propagation results 
(moreover quantitative), which answers partly the question raised by Kac of relating the 
long-time behavior of the A-particle system with the one of its mean-field limit. 
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Finally the general method that we provide for solving these problems is, we hope, 
interesting by itself for several reasons: (1) it is fully quantitative, (2) it is highly flexible 
in terms of the functional spaces used in the proof, (3) it requires a minimal amount of 
informations on the iV-particles systems but more stability information on the limiting 
PDE (we intentionally presented the assumption as for the proof of the convergence of a 
numerical scheme, which was our "methodological model"), (4) the "differential stability" 
conditions that are required on the limiting PDE seem (to our knowledge) new, at least 
at the level of Boltzmann or more generally transport equations. 

1.4. Main results. Without waiting for the full abstract framework, let us give a slightly 
fuzzy version of our main results, gathered in a single theorem. The full definitions of 
all the objects considered shall be given in the forthcoming sections, together with fully 
rigorous statements (see Theorems 13.271 14.11 and 15. ip . 

Theorem 1.1. Let S^ L denotes the semigroup of the spatially homogeneous Boltzmann 
equation (acting on P(M. d )) and , N > 1 denotes the semigroup of the N -particle 
system satisfying a collision Markov process. Then for any (one particle) initial datum 
fo E .P(M d ) with compact support and the corresponding tensorized N particle initial data 
f® N , we have 

(i) In the case where and correspond to the hard spheres collision kernels 

(HS) ; we have for any £ E N* and any N > 2£: 



sup 



sup 



tS(0,oc) ?€ Wrl.»^)W, IMI^l.ocpgd)®* 



(ii) In the case where S?° and correspond to the true (or cutoff) Maxwell 
molecules (tMM)-(GMM), we have for any I E W, any N > 21 and any 

5£(0,2/d): 

Cs 



sup sup 

te(o,oo) ^eJ^.lMI- 



_ (n e [s?(f® N )]-sr(for,<p)<£ 2 

I <1 



2 

Nd- 



for some constant C v E (0, oo) which may blow up when 77 — >• and where 



;i + iei 4 ) \m\ ^< 00 



(p standing for the Fourier transform of ip. 
(hi) In the case where and correspond to the cutoff Maxwell molecules 
(GMM), we have for any £ E N*, any N > 2£, any s E (d/2, d/2 + 1) and any 
TE(0,oo): 



sup 

te(o,T) ^eff 8 



sup 

l A\v\\ H s { 



,<i 



S-, 



N 



(fr)}-s? L (for,^)<t 



s,T 



for some constant C S) t E (0, 00) which may blow up when s — > d/2 or s — >■ d/2 + 1 
or when T — > 00. 

Remarks 1.2. • To be more precise the constants depend on the the initial datum 

fo through polynomial moments bounds for ( GMM) and (tMM), and exponential 
moment bound for (HS). However some one needs also a bound on the support of 
the energy of the N -particle empirical measure. When the N -particle initial data 
are simply tensor products of the 1-particle initial datum, the simplest sufficient 
condition is the compact support of /q. It could be relaxed at the price of an 
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additional error term. This restriction can also be relaxed by considering N -particle 
initial datum conditioned to the sphere S Nd ~ 1 (\^N) of constant energy as in |22j: 
these extensions are considered in Section^ 

• Note that the two first estimates (i) and (ii) are global in time. This is an im- 
portant qualitative improvement over previous chaos propagation results for Boltz- 
mann collision processes. 

• In the third estimate (in), the rate of convergence 0{1/N 1 / 2 ) is optimal in the case 
of Maxwellian kernels, as predicted by the law of large numbers. 

1.5. Some open questions and extensions. 

• What about larger classes of initial data, say with only and as few as possible 
polynomial moments? 

• What about optimizing the rate? Is the rate iV -1 / 2 always optimal and how to 
obtain in general when no Hilbert structure seems available for the estimates (such 
as for the hard spheres case (HS))? 

• What about more general Boltzmann models? A work is in progress for applying 
the method to inelastic hard spheres and inelastic Maxwell molecules Boltzmann 
equation with thermal or stochastic baths. Another issue is the true (without cut- 
off) Boltzmann equation for hard or soft potential: only assumption (A4) (see 
below) remains to be proved in that case. The question of proving uniform in time 
chaos propagation also remains open in this case. 

• What about a central limit theorem and/or a large deviation result with the help 
of our setting? A natural guess would be that a higher-order "differential stability" 
(see below) is required on the limit system, and this point should be clarified. 

• Finally let us mention that works are in progress to apply our method to Vlasov 
and McKean- Vlasov equations. 

1.6. Plan of the paper. In Section[2]we set the abstract functional framework together 
with the general assumption and in Section [3] we state and prove the abstract Theo- 
rem 13.271 In Section 0] we apply this method to the (true) Maxwell molecules: we show 
how to choose metrics so that the general assumptions can be proved (Theorem 14. ID . In 
Section [5] we apply the method to hard spheres molecules and quantitative chaos propa- 
gation in Theorem 15.11 Finally in Section [6] we address several important extensions and 
applications: we relaxe the assumption of compactly supported initial datum by consider- 
ing iV-particle initial data conditioned on the energy sphere as suggested by Kac, and as a 
consequence of our uniform in time chaos propagation result we study the chaoticity of the 
steady state. Finally we conclude with some remarks and computations on the BBGKY 
hierarchy and its link with our work. 



2. The abstract setting 

In this section we shall state and prove the key abstract result. This will motivate the 
introduction of a general functional framework. 

2.1. The general functional framework of the duality approach. Let us set the 

framework. Here is a diagram which sums up the duality approach (norms and duality 
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brackets shall be precised in Subsections 12,31) : 

observables 



E N /& 



N 



Liouville / Kolmogorov 



P (E 

± syin V 



duality 



P N (E) C P(E) 



Liouville / Kolmogorov 



C b (E N ) 
i \ 



Ft 



P(P(E)) 



duality 



W 

Cb(P(E)) 



In this diagram: 

• E denotes a Polish space. 

• (5 denotes the iV-permutation group. 

• P S ym(E N ) denotes the set of symmetric probabilities on E N : For a given permu- 
tation a G & , a vector V = (v\,...,vn) G ^ , a function G C b {E N ) and a 
probability p N G P(E N ) we successively define V CT = (^Vx), Vo-(v)) £ ^ > ^ct 6 
C b (E N ) by setting = (j>(V a ) and G by setting {p% , 0) = (p N , <f> a ). 
We say that a probability on 2? is symmetric or invariant under permutations 



if = p N for any permutation a G 6 



JV 



V 



For any V G £7 the probability measure py denotes the empirical measure: 



N 



PV 



1 N 



V = («i 3 ...,v N ) 



i=i 



where 5^ denotes the Dirac mass on E at point Vj. 

• PjvCE) denotes the subset {/^, V G P^} of P(E). 

• P(P(E)) denotes the set of probabilities on the polish space P(E). 

• Cb (P(E)) denotes the space of continuous and bounded functions on P(E), the 
latter space being endowed with the weak or strong topologies (see Subsection l2.3|) . 

• The arrow pointing from E N /& N to Pn{E) denotes the map defined by 

VV€E N /e N , ^{V):=p^. 

• The arrow pointing from C b {P{E)) to C b {E N ) denotes the following map 7T^ 

V*GC 6 (P(E)), VVeE N , (tt#*) (V) := $ {$) . 

• The counter arrow pointing from C b {E N ) to C\,{P{E)) denotes the transformation 
R N defined by: 



V<l>€C b (E«),Vp€P(E), R N \<P](P) 



In the sequel we shall sometimes use the shorthand notation instead of R\4>\ 
for any I G N* and (f> G C b (E e ). 
• The arrow pointing from P sym (E N ) to P(P(E)) denotes the following transfor- 
mation: consider a symmetric probability p N G P sym (E N ) on E N , we define the 
probability on probability TTpp N G P(P(E)) by setting 

V G G C b (P(E)), = <A*#*> = (P N ,* (Pv)) , 

where the first bracket means (•, -)p(p(E)),c b (P(E)) an d the second bracket means 

(") ~)p{E N ),C b (E N )- 
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• The arrows pointing from the first column to the second one denote the procedure 
of the writing either of the Liouville transport equation associated with the set of 
ODEs of a particle system (as for mean-field limits for Vlasov equations), or the 
writing of the Kolmogorov equation associated with a stochastic Markov process 
of a particles system (e.g. jump or diffusion processes). 

• Finally the dual spaces of the spaces of probabilities on the phase space can be 
interpreted as the spaces of observables on the original systems. We shall discuss 
this point later. 

Remark 2.1. Consider a random variable V on E with law p N G P(E N ). Then it 
is often denoted by pty the random variable on P(E) with law itpp N € P{P{E)). Our 
notation is slightly less compact and intuitive, but at the same time more accurate. 

Remark 2.2. Our functional framework shall be applied to weighted probability spaces 
rather than directly in P{E). More precisely, for a given weight function m : E —> R + we 
shall use (subsets) of the weighted probability space 

(2.1) {p € P(E); M m (p) := (p, m) < oo} 

as our core functional space. Typical examples are m(v) := rh(distE(v, Vq)) } for some fixed 
vq € E and fh(z) = z k or m(z) = e azk , a,k > 0. We shall sometimes abuse notation by 
writing for M m when rh(z) = z k in the above example. 

2.2. The evolution semigroups. Let us introduce the mathematical objects living in 
these spaces, for any N > 1. 

Step 1. Consider a process (Vj ) on E N which describes the trajectories of the particles 
(Lagrangian viewpoint). The evolution can correspond to stochastic ODEs (Markov pro- 
cess), or deterministic ODEs (deterministic Hamiltonian flow). We make the fundamental 
assumption that this flow commutes with permutations: for any a G & N , the solution at 
time t starting from {Vq) is (V/^) . This reflects mathematically the fact that particles 
are indistinguishable. 

Step 2. One naturally derives from this flow on E N a corresponding semigroup acting 
on P sym (E N ) for the presence density of particles in the phase space E N . This corresponds 
to a linear evolution equation 

(2.2) d t f N = A N f N , f N e P sym (E N ), 

which can be interpreted as the forward Kolmogorov equation on the law in case of a 
Markov process at the particle level, or the Liouville equation on the probability density 
in case of an Hamiltonian process at the particle level. As a consequence of the previous 
assumption that the flow (Vj ) commutes with permutation, we have that acts on 
P sym (E N ). In other words, if the law f$ of belongs to P sym (E N ), then for any times 
the law f t N of V t N also belongs to P sym (E N ). 

Step 3. We then define the dual semigroup Tf of acting on functions 4> G C),(E N ) by 

V/ ff eP(B ff ), <f>eC b (E N ), (f N ,T t N (cP)):=(S?(f N )^) 

and we denote its generator by G N , which corresponds to the following linear evolution 
equation: 

(2.3) dt<f> = G N (<t>), cP€C b (E N ). 

This is the semigroup of the observables on the evolution system (V^) on E N . 
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Let us state precisely our assumptions on the iV-particles dynamics. Here and below, for 
a given 1-particle weight function m : E — )■ M + , we define the iV-particle weight function 

1 N 

(2.4) VV = (v 1 ,...,v N )€E N , M m »:=-^raW = ^,m) = M m (^). 

i=l 

Again, we shall sometimes abuse notation by writing Mj? and instead of and 
with m(z) = z k in the example of Remark 12.21 In particular, when E = M. d endowed with 
the euclidian structure, we have 

N 

X 



1 N 

(2.5) VF = K -,vn) e M diV , Mf (V) := - £ N fc = | • | fc > = M fe (p 



i=l 



(Al) On the iV-particle system. Together with the fact that G N and T t N 
are well defined and satisfy the symmetry condition introduced in Step 2 
above, we assume that the following moments conditions hold: 

(i) Energy bounds: There exists a weight function m e and a constant 
£ G (0, oo) such that 

(2.6) supp/f C E N := {V G M^(V) < £} . 

(ii) Integral moment bound: There exists a weight function mi, a time 
T G (0, oo] and a constant C^ mi G (0, oo), possibly depending on T, 
and mi, m e and £ , but not on the number of particles N, such that 

(2-7) sup (f t N ,M^)<C^ mi . 

0<t<T 

(iii) Support moment bound at initial time: There exists a weight func- 
tion 7713 and a constant G (0, +oo), possibly depending on the 
number of particles N and on £, such that 

(2.8) supp/o^ C {V G S^;AC(V0 < Cg^} . 



Remark 2.3. Note that these assumptions on the N -particle system are very weak, and 
do not require any precise knowledge of the N -particle dynamics. 

Step 4- Consider a (possibly nonlinear) semigroup S^ L acting on P{E) associated with 
the limit (possibly nonlinear) kinetic equation: for any p G P(E), S^ L (p) := ft where 
ft G C(R+, P(E)) is the solution to 

(2.9) d t ft = Q(ft), fo = P- 



Step 5. Then we consider its pushforward semigroup T£° acting on Cb(P(E)) defined by: 

vp g p(e), $ g c b (p(E)), rrm(p) ■■= $ {s? L ( P )) . 

Note carefully that Tj°° is always linear as a function of (although of course be careful 
that T t °°[<l>](p) is not linear as a function of p). We denote its generator by G°° , which 
corresponds to the following linear evolution equation on Cb(P(E)): 



(2.10) 



d t <f> = G°°($). 
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Remark 2.4. The semigroup T£° can be interpreted physically as the semigroup of the 
evolution of observables of the nonlinear equation (12. 9p ." 

Given a nonlinear ODE V' = F(V) on M. d , one can define (at least formally) the linear 
Liouville transport PDE 

d t p + V v -(Fp) = 0, 

where p = pt(v) is a time- dependent probability density. When the trajectories (Vt(v)) of 
the ODE are properly defined, the solution of the associated transport equation is given 
by pt(v) = V* t (po), that is the pullback of the initial measure po by V-t (for smooth 
functions, this amounts to pt(v) = po(V-t(v))). Now, instead of the Liouville viewpoint, 
one can adopt the viewpoint of observables, that is functions depending on the position of 
the system in the phase space (e.g. energy, momentum, etc.) For some observable function 
<£>() defined on W^ 1 , the evolution of the value of this observable along the trajectory is given 
by <f>t{v) = 4>o(Vt(v)). In other words we have <fr t = V*(fro = (fro o V t . Then cfr t is solution to 
the following dual linear PDE 

dtcfr-F- V v <fr = 0, 

and it satisfies 

(<fa,Po) = (Vt<f>o,Po) = {(fro,Vl t p ) = {(fro,pt)- 

Now let us consider a nonlinear evolution system V = Q(V) in an asbtract space H. 
Then (keeping in mind the analogy with ODE '/PDE above) we see that one can formally 
naturally define two linear evolution systems on the larger functional spaces P(H) and 
Cb(H): first the transport equation (trajectories level) 

dtir + V- (Q(v) tt) = 0, vr G P(H) 

(where the second term of this equation has to be properly defined) and second the dual 
equation for the evolution of observables 

- Q{v) • V$ = 0, <3? G C b (H). 

Taking H = P(E), this provides an intuition for our functional construction, and also 
for the formula of the generator G°° below ( compare the previous equations with formula 
([2.17p ) and the need for developping a differential calculus in P(E). Be careful that when 
H = P(E), the word "trajectories" refers to trajectories in the space of probabilites P(E) 
(i.e., solutions to the nonlinear equation (j2.9|) ). and not trajectories of a particle in E. 

Another important point to notice is that for a dissipative equation at the level of H, 
one cannot define reverse "characteristics" (the flow in H is not defined backwards), and 
therefore for the nonlinear Boltzmann equation, only the equation for the observables can 
be defined in terms of trajectories. 
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Summing up we obtain the following picture for the semigroups: 

observables 



Pf on E N /& 



If on C b {E N ) 



i4 



P N (E) C P{E) 



observables 





f\ 








R N 




\ 1 




on 


C{P{E)) 



observables 



Sf L on P(E) 

Hence we see that the key point of our construction is that, through the evolution of 
observables one can "interface" the two evolution systems (the nonlinear limit equation 
and the iV-particles system) via the applications tt^ and R N . From now on we shall denote 



7T 



N 



TT, 



2.3. The metric issue. P(E) is our fundamental space, where we shall compare (through 
their observables) the marginals of the iV-particle density / t and the marginals of the 
chaotic oo-particle dynamic f®°°. Let us make precise the topological and metric structure 
used on P(E). At the topological level there are two canonical choices (which determine 
two different sets C(P(E))): (1) the strong topology (associated to the total variation 
norm that we denote by || • ||m0 an d (2) the weak topology (i.e., the trace on P(E) of 
the weak topology M 1 (E), the space of Radon measures on E with finite mass, induced 
by C b (E)). 

The set Cb(P(E)) depends on the choice of the topology on P(E). In the sequel, we 
will denote Cb(P(E),w) the space of continuous and bounded functions on P(E) endowed 
with the weak topology, and Cf,(P(E),TV) the similar space on P{E) endowed with the 
total variation norm. It is clear that Cb{P{E),w) C Cb(P(E),TV). 

The supremum norm ||3>||l°°(p(.e)) does not depend on the choice of topology on P(E), 
and induces a Banach topology on the space Cb(P(E)). The transformations 71"^ and R N 
satisfy: 

(2.11) ||vr 7V $|| i . 00(i?JV) < ||$ || L oo (P(E)) and \\R N [(t>}\\L°°(P(E)) < U\\l™(E")- 

The transformation TT^ is well defined from Cb{P{E),w) to Cb{E N ), but in general, 
it does not map C b (P(E),TV) into C b {E N ) since V £ E N H- $ £ (P(E),TV) is not 
continuous. 

In the other way round, the transformation R N is well defined from Cb{E N ) to Cb(P(E),w), 
and therefore also from Cb{E N ) to Cb{P{E),TV): for any <fi G Cb(E N ) and for any se- 
quence f k -± f weakly, we have f* N f® N weakly, and then R N [(j)\{f k ) -»• R N [<j>]{f). 

The different possible metric structures inducing the weak topology are not seen at the 
level of Cb(P(E),w). However any Holder or C fc -like space will strongly depend on this 
choice, as we shall see. 
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Definition 2.5. For a given weight functions mg : E — > M + and some constant £ G (0, oo), 
we define the subspaces of probabilities: 

(2.12) Pg := {/ G P(E); (f, mg) < oo, (/, m e ) < £), 

where m e is introduced in (Al). For a given constraint function mg : E — > M. D such that 
the components mg are controlled by mg, we also define the corresponding constrained 
subsets 

Pg,r ■■= {f £ Pg; (f,™g) = r}, rGR D , 
the corresponding bounded subsets for a > 

BPg,a := {/ G PgG); (f,mg) < a}, Pg^ a := {/ G BPg^ (f,m g ) = r}, 

and the corresponding vectorial space of "increments" 

TPg := {h -/i;3rGK° s.t. f h f 2 G Pg, r ] 

and IPg, a as expected. Now, we shall encounter two situations: 

• Either distg denotes a distance defined on the whole space Pg{E), and thus on Pg Y 
for any r G M D . 

• Either there is a vectorial space Q D IPg endowed with a norm \\ ■ \\g such that we 
can define a distance distg on Pg Y for any r G MP by setting 

V/, g£Pg,r, distg (/, g) := \\g - f\\g. 

Finally, we say that two metrics do and d\ on Pg are topologically uniformly equivalent 
on bounded sets if there exists k G (0, oo) and for any a G (0, oo) there exists C a G (0, oo) 
such that 

Vf,g€ BPg A d Q (f,g) < C a [dx(f,g)] K , di(f,g) < C a [d (f,g)] K . 

If do and d\ are resulting from some normed spaces G$ and Q\, we abusively say that Qq 
and Q\ are topologically uniformly equivalent (on bounded sets). 

Example 2.6. The choice mg := 1, mg := 0, || • \\g := \\ ■ ||jv/i recovers Pg(E) = 
P(E). More generally on can choose mg k {v) := distE{v,Vo) k , mg k := 0, || • \\g k := 
|| • distE{v,Vo) k \\ M i . For k\ > k2,k^ > 0, the spaces Pg k2 and Pg kz are topologically 
uniformly equivalent on bounded sets of Pg ki ■ 

Example 2.7. There are many distances on P(E) which induce the weak topology, see for 
instance [37J . In section \2. 7| below, we will present some of them which have a practical 



interest for us, and which are all topologically uniformly equivalent on "bounded sets" of 
P{E), when the bounded sets are defined thanks to a convenient (strong enough) weight 
function. 

2.4. Differential calculus for functions of probability measures. We start with a 
purely metric definition in the case of usual Holder regularity. 

Definition 2.8. For some metric spaces Q\ and Qi, some weight function A : Q\ i— )■ 
and some n G (0,1], we denote by C^ v (Qi,Q2) the weighted space of functions from Q\ 
to Q2 with n -Holder regularity, that is the functions S : Q\ — > Q2 such that there exists a 
constant C > so that 

(2.13) V/i./sdi distg 2 (5(/ 1 ),«S(/ 2 ))<C'A(/ 1 ,/ 2 )dist ei (/ 1 ,/ 2 f, 

with A(/i, / 2 ) := max{A(/i), A(/ 2 )} and distg fe denotes the metric of (the tilde sign in 
the notation of the distance has been removed in order to present unified notation with the 
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next definition). We define the semi-norm ['\c a ' r) (g 1 g 2 ) * n ^A^^i'^ 2 ) as ^ e ^ n fi mum °f 
the constants C > such that (|2.13p holds. 

Second we define a first order differential calculus, for which we require a norm structure 
on the functional spaces. 

Definition 2.9. For some normed spaces Q\ and Q2, some metric sets Q\ and Q2 such 
that Qi — Qi C Qi, some weight function A : Q\ 1— > [1, 00) and some r\ G (0, 1], we denote 
C^(Qi,Qi;Q 2 ,Q 2 ) (or simply C^ v (Qi; Q2)) , the space of continuously differentiate func- 
tions from Q\ to Q2, whose derivative satisfies some weighted n-Holder regularity. More 
explicitely, these are the continuous functions S : Q\ — > Q2 such that there exists a contin- 
uous fonction DS : Q\ — > B{Q\,Q2) (where B{Q\,Q2) denotes the space of bounded linear 
applications from Q\ to Q2 endowed with the usual norm operator), and some constants 
Ci > 0, i = 1, 2, 3, so that for any f±, f 2 G Q\: 

(2-14) \\S(f2)-S(f 1 )\\ g2 < C 1 A(f 1 ,f2)\\f2-fi\\g 1 

(2-15) ll(^[/i],/2-/i)|| 0a < C 2 A(f 1 ,f 2 )\\f2-f 1 \\g 1 

(2.16) || t S(/ 2 )_ 5 (/ 1 )_( jDt S[/ 1 ],/ 2 _/ 1 )|| g2 < C 3 A(/!,/ 3 ) ||/ 2 - /iHj^. 

Notations 2.10. For S € C/*, we define Cf , i = 1, 2, 3, as the infimum of the constants 
Q > such that ^T\M) (resp. (pJ5|) . flgTrgp ) holds. We then denote 

[5] c o,i := Cf, [S) c ho := Cf, [S] c i, v := Cf , ||5|| c i,„ := Cf + Cf + Cf , 

and we wiZ remove the subscript A when A = 1. 

Remark 2.11. In £/ie sequel, we shall apply this differential calculus with some suitable 
subspaces Qi C P(E). This choice of subspaces is crucial in order to make rigorous the 
intuition of Griinbaum [20] (see the — unjustified — expansion of Hf in [20\). It is 
worth emphasizing that our differential calculus is based on the idea of considering P{E) 
(or subsets of P(E)) as "plunged sub-manifolds" of some larger normed spaces Qi. Our 
approach thus differs from the approach of P.-L. Lions recently developed in his course 
at College de France [26] or the one developed by L. Ambrosio et al in order to deal with 
gradient flows in probability measures spaces, see for instance [2] . In the sequel we develop 
a differential calculus in probability measures spaces into a simple and robust framework, 
well suited to deal with the different objects we have to manipulate (l-particle semigroup, 
polynomial, generators, ...). And the main innovation from our work is the use of this 
differential calculus to state some subtle "differential" stability conditions on the limiting 
semigroup. Roughly speaking the latter estimates measure how this limiting semigroup 
handles fluctuations departing from chaoticity. They are the corner stone of our analysis. 
Surprisingly they seem new, at least for Boltzmann type equations. 

This differential calculus behaves well for composition in the sense that for any given 
U G and V G C$(g 2 ,G 3 ) there holds S := V o U G C]^{Q U Q Z ) for some 

appropriate weight function A5. We conclude the section by stating a precise result well 
adapted to our applications. The proof is straightforward by writing and compounding 
the expansions of IA and V provided by Definition 12.91 and we then skip it. 

Lemma 2.12. For any given U G C\ r \Q x ,Q2) and V G C l ^{Q 2 ,Qz) there holds S := 
V oU G C^iQuQs) and DS[f] = DV[U(f)\ o DU[f]. More precisely, there holds 

[SU,x < [V] c o,i [UUa, [SU.0 < [V] c i.o [WUo 
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and 

[S] c i,v < [V] C i,o [UU, n + [V] C i,n [uHl 

A 1 +v A o A 

When further V E C l,1 (Q2, G3), we also have S := V o U E Cj^[Qi, G3) with 
[S} c i,n < [V] C i,o [WU, + [V] c i,i [W]J, , (1+J))/2 . 

A 2 A °A 

2.5. The pushforward generator. As a first example of application of this differential 
calculus, let us compute the generator of the pushforward limiting semigroup. Assume 
the following: 

(A2) Existence of the generator of the pushforward semigroup. For 

some normed space Q\ and some probability space Pg 1 (E) (defined as 
above) associated to a weight function nil and constraint function mi, 
and endowed with the metric induced from Qi, for some 5 E (0,1] and 
some a E (0, 00) we have for any a E (a, 00): 

(i) For any t E [0, +00), Sj^ L : BPg lya — > BPg x ^ a is continuous (for the 
metric distc^), uniformly in time. 

(ii) The application Q is 5- Holder continuous from BPg lj0 into Q\. 

(hi) For any / E BPg 1}CL , for some r > the application [0, t) — > BPg 1}CL , 
t S t NL (f) iS C^^T^BPg^g), with S(f)'(0) = Q(f). 

Lemma 2.13. Under assumption (A2) the generator G°° of the pushforward semigroup 
Tj°° exists as an unbounded linear operator on C(BPg 1:a (E);M) for any a E (a, 00) with 
domain including C 1 ' S (BPg 1:a (E);M.), and on this domain it is defined by the formula 

(2.17) V$6^(BP Slia (£);l), V/ E BPg 1>a (E), (G°° $)(/):= (!>*[/], Q(/)) . 

Proof of Lemma \2.13[ We split the proof in several steps. 

Step 1. First (Pf ) is a continuous semigroup on C(23Pg lia (P);R): consider $ E C(BPg ua {E) 
and some sequence (/„) of BPg 1>a (E) such that distg 1 (/ n , /) — > 0, then thanks to (A2)-(i) 
we deduce (7f $)(/„) = $(Sf L (jW)) -> $(Sf L (/)) = (P t °°$)(/), and by composition 
2f$€C(BP 0l , a (E);R). 
Next, we have 

= sup II = sup sup \$(S t NL (f))\<l, ||*||= sup |$( 5 )|. 

II*II<1 ||$||<i /eBP gi , a seSPg lia 

Now, from (A2)-(iii) there exists a modulus of continuity w such that \\S^f — f\\g 1 < 
uj{t) — > as t — > for any / E £>Pg li(1 (e.g. ui(t) = t a , a E (0, 1)), which implies 

V*€C(BP 0l , a (S);R), ||rr*-$||= sup |<D(Sf L (/)) "$(7)1^0. 

/6BP Sl ,a 

We therefore deduce from Hille-Yosida's Theorem that (Tf°) has a closed generator G°° 
with dense domain included in C(HPg lj0 (P);R). 

Step 2. Let us define G°°$ by 

V* E C^OBPc^R), V/ E BPg ua , (G°°$)(/) := <£>$[/], Q(/)) . 

The right-hand side quantity is well defined since D$>(f) E S(£7i,R) = and Q(/) E 
C/i. Moreover, since both / 1— >• D$>[f] and / 1— > Q(f) are continuous we have G°°$ E 
C(BPg 1)0 ;R). 
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Step 3. Consider G C 1,<5 (i3Pg 1>a ; R). By the composition rule of Lemma 12.121 for any 
fixed / G P 5l (E), t r t °°$(/) = $ o S? L (f) is C 1 ' <5 ([0, r/); M) and 

|er*)(/)U = s(*°*(/)<'))L, 

W(/))U4 5(/) I« 

= (£)*[/], Q(/)) = (G°°*)(/), 
which precisely means that any such $ belongs to domain(G°°) and that (12. 171) holds. □ 

2.6. Compatibility of the ir N and Rr transformations. Our transformations tt n and 
Rr behave nicely for the sup norm on Cb(P(E),TV), see (|2.1ip . More generally we shall 
consider "duality pairs" of metric spaces: 

Definition 2.14. We say that a pair {J-,Q) of normed vectorial spaces are "in duality" if 

(2-18) V/ £ g,V^J < WfWgWWr 

where \\ ■ \\jr denotes the norm on T and \\ ■ \\g denotes the norm on Q. 

The "compatibility" of the transformation R? for any such pair follows from the mul- 
tilinearity: if J- and Q are in duality, T C Cb(E) and Pg is endowed with the metric 
associated to || • \\g, then for any ip = (p\ X • • • X (pi E J-®^, the polynomial function Rr is 
of class C 1,,? (Pg,IR) for any n G (0, 1]. Indeed, defining for f a E Pg 1 , a = 1,2, 



g Rj hh> := E II / ^ 

i=i v^-^ / 

we have 

<(/ 2 )-^(/i) = e( n / (w,/ 2 -/i) ( n / ^rf/i), 

and 

4(/ 2 ) - 4(A) - £><[/i](/ 2 - /i) = 

= e ( n / ) {<Pjj2-h) ( n / ^dfi) (vi,/ a -/i)( n / 

We deduce then for instance G C 1 ' 1 (^;M) since 

\R e v (f 2 ) - < * IbH^^),-! ||/ 2 - /lib, |^<[/l] (h)\ < i M\ T ^y-r \\h\\g 

(2.19) |4(/ 2 ) - 4(/0 - DR^[h]{f 2 - h)\ < £ -^-H\\ ip \\^ {LOO)e , 2 ||/ 2 _ ^115, 
where we have defined 

IMI.F*®CL*>)'-* := . . max ri1 „ n \\<Pii\\r ■ ■ ■ \Wi k \\T I [ \Wj\\h°°(E)- 

distincts in [|l,t|J 

Remarks 2.15. • It is easily seen in this computation that the tensorial structure 

of (p is not necessary. In fact it is likely that this assumption could be relaxed all 
along our proof. We do not pursue this line of research. 
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• The assumption T C C b (E) could also be relaxed. For instance, when T := 
Lipo(-E) is the space of Lipschitz function which vanishes in some fixed point 
xo G E, Q is its dual space, and Pg := {/ G P\{E); (/, diste(-, Xq)) < a} for 
some fixed a > 0, we have R G C 1,1 (Pi(E); R) wf/i 

ro^i ^ e„£-i ii, oil rp^i 1(1 — 1) e_i ■■ ■■ 

[-K^JcO' 1 < « 1Mb®*, [-K^Jai.i < — 2 — llVlb®^ 

or equivalently R*, G C^ ,:L (Pi(-E); R) ictifc A(/) := . 

In the other way round, for the projection tt n it is clear that if the empirical measure 
V i — y Hy belongs to C k,v (E N ,Q) for some norm space Q G M l (E), then by composition 
one has 

(2-20) ll 7riY ( $ )!lc fc >"(s JV ;K) - ^ H $ llc fc ."(6) • 

However the regularity of the empirical measure depends on the metric Q. 

Example 2.16. In the case T = (C b (E),L°°) and Q = (M 1 (E),TV), (bT20j) is trivial 
with k = r) = 0. 

Example 2.17. When T = Lip (E) (Lipschitz function vanishing at some given point 
vq) endowed with the norm \\4>\\Lip and Pg{E) (constructed in Example \2.19]) is endowed 
with the Wasserstein distance W\ with linear cost, one has fl2.20f) with k = 0, rj = 1: 

-3 fay) | ^ Pllco,i(P 6) Wi(/4,£y) < \m\ c °,HP g )\\X-Y\\ e i, 
where we use A2.22\) . which proves that 

ll 7r7V (^)llc°. 1 (S JV ) - II $ II C^iPg), 

when E N is endowed with the £ l distance defined in l2~i 



2.7. Examples of distances on measures. Let us list some well-known distances on 
P(R d ) (or on subsets of P(W 1 )) useful for the sequel. These distances are all topologically 
equivalent to the weak topology a(P(E),Cb(E)) (on the sets BPk, a {E) for k large enough 
and for any a G (0, oo)) and they are all uniformly topologically equivalent (see [1HE] an d 
section 12.81) . 



Example 2.18 (Dual-Holder (or Zolotarev's) distances). Denote by distE a distance on 
E and let us fix vq G E (e.g. vq = when E = R rf in the sequel). Denote by C ' S (E), 
s G (0,1) (resp. Lip (E)) the set of s -Holder functions (resp. Lipschitz functions) on E 
vanishing at one arbitrary point v$ G E endowed with the norm 

[<f]s := sup l^)-^)l ;gg (0 !] [ ] := Ml . 
X ^ E <±\st E {x,y) s 

We then define the dual norm: take mg := 1, mg := and Pg(E) endowed with 

(9 ~ f, f) 



(2.21) V/, 9 e? 5 , [g-f]*-= sup 



[<p]t 



Example 2.19 (Wasserstein distances). For q G (0,oo), define W q on 
Pg(E) = P q (E) := {/ G P(E); (/, dist(;V ) q ) < oo} 

by 

Vf,geP q (E), W q (f,g):= inf / dtst E (x, y y U(dx,dy), 

nen(/, 5 ) J E xE 
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where Il(f,g) denote the set of probability measures II G P(E x E) with marginals f and 
g (n(A,E) = f(A), U(E,A) = g{A) for any Borel set A C E). Note that for V,Y £ E N 
and any q G [l,oo), one has 



(2.22) W q (i4,fj#) = oI £He n /6n) (V,Y) := mm (± £ dist E (v u y a{€) )^ 



1/9 



and that 

(2.23) VfigeP^E), W 1 (f,g) = [f-g}{= sup (f-g,<P). 

4>£Li Po {E) 

We refer to [46J and the references therein for more details on the Wasserstein distances. 
Example 2.20 (Fourier-based norms). For E = M. d , mg 1 := \v\, mg 1 := 0, let us define 

V/€TP 6l , ||/|| 6l = 1/1. := sup ffi, s€(0,l]. 

^eM d Is I 

Similarly, for E = M. d , mg 2 := \v\ 2 , mg 2 := v, we define 

V/GTPg 2 , H/llfe = |/|. := sup s€(l,2]. 



Example 2.21 (More Fourier-based norms). More generally, for E = M. d , k G N* ; we 
define rag := \v\ k , m g := (v j ) jm d t \j\< k -i where for j = (ji,...,jd) G N d we set v j = 
{v[ l ,...,v 3 d d ) and \j\ = ji + ■■■+ 3d, and 

V/GTPg, \\f\\g = \f\ s := sup a € (0,k\. 



In fact, we may extend the above norm to Ml(R d ) in the following way. We first define 
for f G Ml_ 1 (R d ) and j G N d , \j\ <k — l, 

Mj[f]--= [ v*f(dv). 

For a fixed (once for all) function x in the Schwartz space S(R d ) such that x = 1 on the 
set {v G M. d , \v\ < 1}, so that in particular j Rd J-^ 1 (x)(v) dv = x(0) = 1, we define M k [f] 
be its Fourier transform 

M k [f}(!;) := X (0 I £ M Af\^i m 
\\j\<k-l 3 ' 

which is some smooth version of the Taylor expansion of f at £ = 0. Then we may define 
the seminorms 

|/|fc == sup (\i\- k f(0-M k [f}(0) 

and 

:=\f\ k + Y, \ M M- 



j£N d , \j\<k-l 
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Example 2.22 (Negative Sobolev norms). For any s G (d/2,d/2 + 1/2) take E 
mg 1 := \v\, mg 1 := and 



V/gTP 6i , 



9i 



no 



L 2 



Similarly, for any s G [d/2 + 1/2, d/2 + 1) take E = W 1 , mg 2 := \v\ 2 , mg 2 := v and 



Vf£TPg 2 



no 



L 2 



2.8. Comparison of distances when E = M. d . 
Lemma 2.23. Let f,g G P(M d ), £/ien 

(2.24) Vg,fe6 (l,oo) #) < W q (f,g) < M\~^ W\(f , g) c 

with a := 1 — (q — l)/k, 

V s G (0,1], < < (/,<?), 



(2.25) 

(2.26) 

(2.27) 
with 

(2.28) 
with 

and 



V«€(d/2,d/2 + l), <C|/- 5 |? s - rf , C7 = C(d, S )>0, 

Vs>0, A:>0 [/-^^CM^lZ-ffU 1 , C = C(d, s, fc) > 0, 

d 



71 == 



d + fc + /c(d + s-l)' 11 ' d + k + k(d + s-l) 
Vs>l,fc>0, [/- < CM^ ||/ -5[|J_„ C7 = C7(d, S ,A:) > 0, 

d/2 A; 



a 2 : = 



d/2 + k + k(s — I)' 



72 - 



d/2 + k + k(s — 1) 



M k :=max\ f (1 + \x\ k ) f(dx) ; [ (1 + \x\ k ) g(dx) ) . 

Proof of Lemma \2.23\ We split the proof in several steps. For the proof of (|2.24p we refer 
to pi El. 

Proof of (12351) . Let vr G U(f,g). We write 



1/(0-3(01 



< 



< a 



- iv < -e~ iw <)ir{dv,dw) 
lv < -e- iw -t\ir(dv,dw) 
\v - w\ s n(dv,dw), 



which yields (|2.25p by taking the supremum in £ G W d and the infimum in tt G n(/, 5). 
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Proof of (pT26l) . Consider R > and the ball B R = {x G R d ; \x\ < R}, and write 



\\f-9\ 



H~ 



< \f~9\lf 
Je 



\a\ 2s 

Br lei 2 ^ 



Bi 



+ 4 



\t\ 2s 



2,s 



R 

,d-2s 



< C(d)R d - 2(s - 1 '>\f -g\l + 4R a 
Then (12,261) follows by choosing (the optimal) R := \ f — g^ 1 . 

Proof of (I2.27p . We introduce a truncation function xr( x ) = x( x /R), R > 0, where 
X G C°°(R d ), [x\i < 1, < x < h X = 1 on 5(0,1), suppx C 5(0,2), and a mollifer 
function u> £ (x) = e~ d oj(x/e), e > where for instance w(x) = (27r) _d//2 exp(— |x| 2 /2) (and 
thus u e (£) = w(e£) = exp(-e 2 |£| 2 /2)). Fix <p G W 1 ' 00 ^) such that [v?]i < 1, p(0) = 0, 
define tp R :=tpxR, VR,e = Vr*u £ and write 



V (df ~ dg) 



PR,e (df ~dg)+ I (ip R - <p R}E ) (df - dg)+ \ (ip- cp R ) (df - dg). 



For the last term, we have 



(2.29) 



(<PR ~ <P) (df ~ dg) 



< 



< 



1(1- X r) V (df + dg) 



Mi 



ifc+i 



Bi 



R k 



(df + dg) < 



M k+ i[f + g] 
R k 



where Mk + i[f + g] denotes the (k + l)-th moment of / + g. In order to deal with the 
second term, we observe that 

|V^| < X (x/R) + M |V(xh)| < x(x/R) + ^ \V X \(x/R), 

so that for any q G [l,oo] there holds || V^r\\li < CR d / q , for some constant depending 
only on X) d. Next, using that 



we find 
(2.30) 



Wr - <PR,s\\oo < ||V^ii||oo / uj £ (x)\x\dx < C e, 



< Ce. 



(<PR ~ <Pr,e) (df - dg) 



For the first term, using Parseval's identity, 

1 

2^ 
1 



¥r,s (f - g) 



< 2^l|V^|| L i 



¥RU £ (f - g)d£ 
f-9 



izr 1 eM-e z \emdt 



2 \i-\2 



< CR a 



1 + \y\)x(y)dy) \f-g\ 8 s 



-(d+s-l) 



is— 1 



e 2 dz 



(2.31) 



< C R d e- {d+s - l) \f 
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[/ - 9\\ <c(e+ Mk+1 ^ + 9] + R d e-^ \f - g^. 
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This yields (|2.27p by optimizing the parameters e and R. 

Step 4. Proof of (|2.28[) . We start with the same decomposition as before: 



ip (df -dg) = J ip R)£ (df -dg) + J (ip R - ip R>e ) (df - dg) + / (cp - ip R ) (df - dg). 
The first term is controled by 

(f-g) 



J 



<PR,e (df - dg) 



with 



1/2 



(2.32) 



The second term and the last term are controled as before by C e and M^ + i R k respec- 
tively. Summing we obtain 

[/ - 9]\ <c(e+ Mk+l g + 9] + B*" e -(-*> |/ - g\ 6 J . 

This yields ()2.28p by optimizing the parameters e and R. □ 

2.9. On the law of large numbers for measures. 

2.9.1. Remark on the meaning of the W function. For a given function D : P(E) x P(E) — > 
R + continuous for the weak topology, and such that D(f,g) = if and only f = g (D 
stands for a distance on P(E) or a function of a distance), we define 

WE(f):=[ D(^-f) f® N (dV). 
Je n 

On the one hand, from the definition of TTpf® N in section [2TTT we have 
<(/) = W# (/»";/) = C /) , 

with 

W£ (/*;/):= / D(^,f)f N (dV) 

J E N 



Vfe P(E N ) 

\/tt€P(P(E)) W%fcf):= [ D(p, f) n(dp). 

JP(E) 
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2.9.2. Comparison of functions for different distances D. Let us begin with an ele- 
mentary result. 

Lemma 2.24. If D\ < C for some constants C > and r < 1, we have for some 
constant C = C'(C, r) and for any f G P(E) 

(2.33) VVft(/)<C (><(/)) r . 

Proof of Lemma \2.24\ For any e > 0, we have from Young's inequality 

<(/) < C / D 2 (^;/) r /^(dF) 

[ ( i_ r)£ r/(i-) + !: jD2 (^ ;/ )i f ®N (dV ) 



< c 



< C{r) 



r/(l—r) + WUf) 



This yields (|2.33p by optimizing e > 0. 

Lemma 2.25. We have the following rates for the W function: 

• For any f G P2(M. d ), any s G (d/2, d/2 + 1) and any N > 1 i/iere ZioZds 



□ 



(2.34) 



(2.35) 



(2.36) 



I II ZT-. 



(/) 



< Cst(d, M^iV- 1 . 



For any n > i/iere exists k > 1 sitc/i i/iai /or any / G P^(]R ) and any iV > 1 
there holds 

^Cst^^M^iV- 1 /^) . 

For any n > i/iere exists k > 2 snc/i t/iai /or any / G -Pfc(R d ) and any TV > 1 
i/iere /io/ds 

W^ 2 (/)<Cst(n,£;,M,0A^ 1/{d+ ^ 



Remark 2.26. Estimate (|2.36|) ftas to oe compared with the following classical estimate 
established in |37j : /or any / G -P,j+5(R d ) and any N > 1 there holds 



(2.37) 



W$ 2 (/) < Cst(d,M (i+5 )iV-^4. 



zs worth mentioning that (|2.36p improve (|2.37p w;/ien d < 3 and A; is Zarye enough (so 
that 7] <2-d/2). 

Proof of Lemma \2.25[ We split the proof into two steps. 
Proof of CLMD. Let us fix / G P 2 QR d ). First, writing 

i ^ 



we have 



(/) 



iV 2 



AT 

£ 

ii j2=i 



d^j /^(dF) 
(e^r* - /(£)) (e~^ -/(£)) 



d£/ 0jv (dF). 



]R(JV+l)<i 



2 s 
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We use 



and 



f(O)f(dv j ) = 0, j = l,...,d, 



-i v-£ 



/(0 f(dv) 



1 _ e -iK _ e --« + f{dv) 

i-l/(Ol a , 



to deduce 



A' 



AT2 
1 

iV 
1 

iV 



0=1 



d£f® N (dV) 



'N+l)d 



f(0 



|£|2. 

Ci-l/COl 2 ) 



2 



lei 



2 s 



Finally, observing that /(£) = 1 + i (/, v) • £ + 0(M 2 |£| 2 ), and therefore 

l/(0i 2 = (i + i(f,v) ■ z + o(m 2 i^i 2 )) (i - i (f,v)- e + o(m 2 iei 2 ; 
= i + o(M 2 iei 2 ), 

we obtain 



1 

N 



M 2 



iei<i 



|£|2 (« 



iei>i 



(i-l/(C)l 2 ) M 
\i\>i l^l 2s , 



from which (|2.34p follows. 

Proof of ([235]) and ([236]) . By gathering flSgD , (l23H|) in Lemma EES] and Lemma ETJl 
we straightforwardly deduce 



<8>iV/ 



< c7 M ,,(M fc+1 ) /" (n^r - /n^ s ) 72/2 

< C k!S>d (M k+1 )N-^ 2 , 

from which we deduce (|2.35f) because in the limit case k = oo we have 72/2 = l/(2s) and 
we may choose s as close from d/2 as we wish. 

Proof of 0736]) . Estimate (pZ736l) follows from Q735]) with the help of Lemma E723] and 
(MP . □ 
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3. The abstract theorem 



3.1. Assumptions for the abstract theorem. Assume that 

• (Al) and (A2) hold, so that in particular the semigroups , , S^ L and T + °° 



are well defined as well as the generators G N and G 



(A3) 



(3.38) 



Convergence of the generators. In the probability metrized set Pg 1 intro- 
duced in (A2) (associated to the weight function mg 1 and constraint function 
mgj we define 

K Gl := {rGM D ; 3 / G Pg x s.t. vag 1 (/) = r}. 

Then for the weight function Ai(/) = (/, mi) and for some function £2(N) 
going to as N goes to infinity, we assume that the generators G N and G°° 
satisfy 

V$G fl C^(P Ql y,R) 

M^y^TTN-^G^) $ 



< e 2 (N) sup [$U,„ (p 

r6R Sl A i ( 61 ' 



where was defined in (12. 4h . 



(A4) Differential stability of the limiting semigroup. We assume that the 
flow Sj is Cj^ n (Pg lir , Pg 2 ) for any r G Hg ± in the sense that there exists 



> such that 



(3.39) 



sup 

reR ei Jo 



Pt J cl;"(P ei , r ,P 62 ) 



, roNL]l+v' 

+ [t J ^f(P ei . r ,P 52 ) 



where 77 G (0, 1) is the same as in (A3), (rf ,rj") = (77, 1) or (rf,rj") = (1, (1 + 
r})/2), A2 = A^ 1 " 1 "'' ^ and Pg 2 = {p, p G -P^} but it is endowed with the norm 
associated to a normed space Q2 D 



(A5) Weak stability of the limiting semigroup. We assume that, for some 
probabilistic space Pg 3 (E) (associated to a weight function mg 3 , a constraint 
function mg 3 and some metric structure distg 3 ) and that for any a,T > 
there exists a concave and continuous function Q a ,T ■ ^+ — > ^+ such that 
@a,r(0) = 0, we have 

yfi,f 2 eBPg 3 , a (E) 

(3.40) sup distg 3 5f L (/ 2 )) < 6 a , T (dist 63 (A, / 2 )) . 

[0,T) 



3.2. Statement of the result. 

Theorem 3.27 (Fluctuation estimate). Consider a family of N -particle initial conditions 



/q G P sym (E ), N > 1, and i/ie associated solution f t 



N 



&t fo ■ Consider a l-particle 



initial condition /q G P(E) and the associated solution ft = fo- Assume that (Al)- 
(A2)-(A3)-(A4)-(A5) hold for some spaces Pg k , Qu andT k , k = 1,2,3 with Ft C Cb(E), 
and where Tk and Qj. are in duality. 

Then there is an explicit absolute constant C G (0, 00) such that for any N, £ G N*, with 
N > 11, and for any 



tp t g (JinJ 2 nJ 3 ) 
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we have 
(3.41) 



sup 

[0,T) 

< c 



(s?(fo N )-(S? L (fo)) m ),v) 



+£ |M 



'>r<N 

U 0,m 3 > 



t ( w dist c i^P./oV/o 



S g 



where W^istg stands for the Monge-Kantorovich distance in P(Pg 3 (E)), see example \2.19l 
which means for such particular probabilities 

(3.42) W dist (tt£/ V /o ) = / die^ 8 (^/o)/o"W- 

Remark 3.28. Our goal here is to treat the N -particles system as a perturbation (in 
a very degenerated sense) of the limiting problem, and to minimize assumptions on the 
many-particle systems in order to avoid complications of many dimensions dynamics. 

2) In the applications the worst decay rate in the right-hand side of (|3.41|) is always the 
last one, which deals with the chaoticity of the initial data. 

3) It is worth mentioning that in the case when f^ = f® N we have 

and the decay rate is obtained thanks to the law of large numbers for measures presented 
in section \2lh For more general initial datum we refer to section® 

3.3. Proof of Theorems 13.271 For a given function ip G (J 7 ! D T% PI T3) , we break up 
the term to be estimated into four parts: 



< 
+ 
+ 



S?(f N ), v ®l 



®N-e 



Lp &> 1 

qN 1 fN 

a t uo 



®N-l 



< 



Jo > T t (^°My) 



/ O w ,(r<)o^; 



(s?>(fo)r e ,<p 



(=■■ 71) 
=:T 8 ) 
(=■■%). 



We deal separately with each part step by step: 

• 71 is controled by a purely combinatorial arguments introduced in [20 . In some 
sense it is the price we have to pay when we use the injection ir N ; 

• 72 is controled thanks to the consistency estimate (A3) on the generators, the 
differential stability assumption (A4) on the limiting semigroup and the moments 
propagation (Al); 

• is controled in terms of the chaoticity of the initial data thanks to the weak 
stability assumption (A5) on the limiting semigroup and (Al)-(iii). 

Step 1: Estimate of the first term 71. Let us prove that for any t > and any N > U 
there holds 

2 ^ 2 IMIl°°(£*) 



(3.43) 



71 := 



.N 



< 



N 



Since (/q ) is a symmetric probability measure, estimate (|3.43p is a direct consequence 
of the following lemma 
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Lemma 3.29. For any ip € Cb(E e ) we have 



(3.44) 



ViV > 2£, 



(tp® I s 



sym 



< 



2 ^ 2 II^IIl°°(^) 



iV 



where for a function <fi £ Cb (E ) , we define its symmetrized version 4> sym as: 

(3.45) 4> sym = j-^—r y~] 4> a . 

Iojv ^ 



^4s a consequence for any symmetric measure f N € P(E N ) we have 

2 ^ 2 IMIl°°(£«) 



(3.46) 



(r,K(^))-(r^) 



< 



Proof of Lemma 13.291 For a given £ < N/2 we introduce 

A N £ := Uii,...,ii) G [|l,iV|]* : VA; ^ A/, i fc ^ i fc / j and jB^ := A^. 

Since there are N(N — 1) . . . (N — I + 1) ways of choosing t distinct indices among [|1, 
we get 



\B 



^)- l -(S-H)) 



< 



1 -«p(- 2 E^N]v. 



i=0 



where we have used 

Vac €[0,1/2], ln(l-x) > -2x and Vi 6 R, e~ x > 1 - x. 
Then we compute 



1 * 



(n,...,t e )eA Nt i 
1 1 



,Vi 



(h,...,ii)eB N>i 

•2 



N e (N - £) 



Hi S ^(^(i),---,^ W ) + ° 



TV 



and the proof of (j3.44p is complete. Next for any f N G P(E N ) we have 

W = (^(»«'"t), 

and (pTlB"]) trivially follows from (f3l3]) . 
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Step 2: Estimate of the second term 72- Let us prove that for any t G [0, T) and 

any N > 21 there holds 



(3.47) 



T-2 



< Qr.ma Ct IMIoo,.F 2 2 ®(£°°)'- 2 ^ e W- 

We start from the following identity 

/V. [ 

From assumptions (Al) and (A3), we have for any t G [0,T) 



Tf n N - 7T N Tr = - I f g (T t N _ s n N T?) ds = I T" K \G N n N - vr^G 00 ] ds. 



< ( sup (^,<)) (7 



i°°(Ejv) 



ds 



(3.48) 



<s(N)Cff m sup 



reR 5l JO 



T 



s K <p 



ds. 



Now, let us fix r G R gi . Since T t °°{R e v ) = flj o Sf L with Sf L G C^(Pg lir ; Pg 2 ) thanks 
to assumption (A4), and G C 1,1 (Pg 2 ;M) because </? G J 7 ®^ (see subsection I2.6p . we 
obtain with the help of Lemma EjJ that T^(R e J G C 1 '? $ ,(P gi y,R) with 



-iNL~l 



JVL1 

2 



R l 



1/ (l-\-0 f ) 

, (|2.19p and assumption (A4), we hence deduce 
(3-49) jf [Tr«)] c i,, (Psi r) d S < Cfl 2 (\M W oo Y -2 + IMI^^-i) 



With the help of A 2 = A{ 



Then we go back to the computation (|3.48p . and plugging f|3.49f) we deduce f|3.47j) . 
Step 3: Estimate of the third term 73. Let us prove that for any t > 0, N > £ 



(3.50) T 3 :~- 



fF,(T t °°Rt)o^ 



.N 



S?(fo)) ,<p 

< [Ripley c c f m3 ,r ( w i,^g 3 (^p/oV/o)) ■ 



We shall proceed as in Step 2, using that: 

• SUpp7Tp fX c JC := {/ G Pg 3 ; M ma (f) < C^ m3 } thanks to assumption (Al), 

• Sj^ L satisfies some Holder like estimate uniformly on K, and [0, T) thanks to as- 
sumption (A5), 

• R[p G C^iPg^R) because G jf. 
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We deduce thanks to the Jensen inequality (for a concave function) 

r 3 = 



fo N ,K (sfW))) - {fo N ,K (^ L (/o)) 

/( f,<(^(^))-<(^(/o)) 



< 


[ R v. 


c° 


"(Pes) 


(fo N 


dis^ (s? L (fo),s? L (»m 


< 


[Rip_ 


c° 




(fo N 


e ajT (distg 3 (/ ,Mv))> 


< 


[R v _ 


c° 




e a , T 


«/ ", distft,(/ ,A#)». 



Now, by definition of the optimal transport Wasserstein distance we have 

V Hi, (i 2 € P (Pg 3 ) , Wi,p g (m,fJ>2) = inf / distg 3 (//,//)7r(d/i,d)u')> 

7ren(/ii,At 2 ) JPg 3 xPg 3 

where H(nx, M2) denotes the probabilty measures on the product space Pg 3 x Pg 3 with first 
and second marginals fJ,x,fJ-2- in the case when 1x2 = Sf then H((j,x,8f ) = {fix &> <5/ } has 
only one element, and therefore 



Wl,Pg 3 (^p/oV/ ) = mf / dist 03 (/,2)7r(d/,d<7) 

7ren(7r^/,f ,5 /o ) Jp^p^ 

di S tg 3 (f,g)TT P ! fS f (dg)5 f0 (df) 

distg 3 (/ ,5') npfo(dg) 
distg 3 (^J )fo N (dV). 



'Pg 3 xPg 3 

uistg 3 <J0,yj 7 

'P S3 

,JV * -\ f N 1 

uisug 3 VA 

We then easily conclude. □ 

4. (True) maxwellian molecules 

4.1. The model. Let us consider E = M. d , d > 2, and a iV-particles system undergoing 
space homogeneous random Boltzmann type collisions according to a collision kernel B = 
T(z) 6(cos 9) (see Subsection II More precisely, given a pre-collisional system of velocity 
variables V = (v\, ...,ujv) £ ^ = (R ) , the stochastic process is: 

(i) for any i' ^ j' , draw a random time T r (|„.,_^.,|) of collision accordingly to an 
exponential law of parameter r(|ity — Vji\), and then choose the collision time Tx 
and the colliding couple (vi,Vj) (which is a.s. well-defined) in such a way that 

Tx = T T(lv .^ Vjl} := ^mm^T r(K ,_ v .,|); 

(ii) then draw a £ S^ 1 according to the law fr(cos Oij), where cos 9™ = a-(vj — Vi) /\vj — 
Vi\; 

(iii) the new state after collision at time T\ becomes 

V*j = (vx,...,v*,...,v*,...,v N ), 

where only velocities labelled i and j have changed, according to the rotation 

-t \ * Vi + vj \vi-vj\a Vi + vj 

(4-1) «i=— 5-^- + 5"^. v j = ^T^ 
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The associated Markov process (Vt) on the velocity variables on (R d ) N is then built by 
iterating the above construction. After scaling the time (changing t —> t/N in order that 
the number of interactions is of order 0(1) on finite time interval, see |39j) we denote 
by the law of Vj, the associated semigroup, G N and respectively the dual 
generator and dual semigroup, as in the previous abstract construction. The so-called 
Master equation on the law is given in dual form by 

(4-2) d t (f t N ,<p) = (f t N ,G N <p) 

with 

i - r 

(4.3) (G N <p)(V) = jj ^r(h-^-l) J d J{ C o S 6 tJ ) [tplj-tp] da 

where = ip(V$) and if = <p(V) G C b (R Nd ). 

This collision process is invariant under velocities permutations and satisfies the micro- 
scopic conservations of momentum and energy at any collision time 

N N N N 

= and in 2 = £Ki 2 = £N 2 = m 2 . 

As a consequence, for any symmetric initial law f§ G P sym (R Ard ) the law at later times 
is also a symmetric probability, and it conserves momentum and energy: 

where (fj, a )i< Q <d denote the components of € and 

(4.4) V0:M+^IR + , / <P(\V\ 2 ) f t N (dV) = [ <P(\V\ 2 ) f» (dV) 

JR dN JR dN 

(equality between possibly infinite non-negative quantities). 

The (expected) limiting nonlinear homogeneous Boltzmann equation is defined by (jl.2p , 
(j 1 . 3 p . (|1.4|) . The equation generates a nonlinear semigroup S^ L (fo) := ft for any /o 6 
P2(^ d ) (probabilities with bounded second moment): for the Maxwell case we refer to 
|421I44|. for the hard spheres case we refer to [35] in a L 1 setting, and for the generalization 
of these L 1 solutions to P2(R d ) we refer to [II], [17] and [28]. For these solutions, one has 
the conservation of momentum and energy 

Vt>0, / vf t (dv)= [ vf (dv), f \v\ 2 f t (dv)= f \v\ 2 f (dv). 

JR d JR d JR d JR d 

Without restriction, by using the change of variable a i— )■ —a, from now on we restrict 
the angular domain to 9 G [—tt/2, 7t/2] for the limiting equation as well as the iV-particle 
system. Therefore we assume supp6 C [0, 1]. We still denote by b the symmetrized version 
of 6 by a slight abuse of notation. 

4.2. Statement of the result. In this section we consider the case of the Maxwell 
molecules kernel. More precisely we shall assume 

(4.5) T = l, b G L,~([0,1)), Vt?>1/2, CJb) := [ b(cos6) \l-cos 9\ v / 2 da < oo. 

J§d-1 
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Indeed for any positive real function ip and any given vector u G M. d we have 



ip(u • a) da 



ip (cos 9) sin 



d-2 



>d6 



's d - x JO 
and the true Maxwell angular collision kernel (tMM) defined in Subsection 11.11 satisfies 
(in dimension d = 3) b(z) ~ K (1 — z)~ 5 / 4 as z — >• 1, which hence fulfills (14. 5p . These as- 
sumptions also trivially include the Grad's cutoff Maxwell molecules (GMM) introduced 
in Subsection ll.il 

Our fluctuations estimate result then states as follows: 

Theorem 4.1 (Boltzmann equation for Grad's cut-off/true Maxwell molecules). Let us 

consider an initial distribution Jq G P(M. d ) with compact support. Let us consider a hi- 
erarchy of N -particle distributions f^ = SP{f^ N ) issued from the tensorized initial data 
Then for any 



fN 

Jo 



JO ■ 



with 

we have 
(4.6) 



sup 

t>0 



< C 



■■= I (i + iei 4 )i^(e)i^<oo), 

Jm d J 
S t N (f N )-(S? L (fo)f N ' 



f a \W\ 



N 



< 



1,2 C V 



+£ \\<p\ 



\f\\jr(git 



w 1 



i(L°°)«- 1 Ww 2 (/o) 



for some constants C^ ^, C v , G (0, oo) which may blow up when r] — > (and depend on b 
and£ :=M 2 (f )). 

In order to prove Theorem 14. 11 we have to establish the assumptions (A1)-(A2)-(A3)- 
(A4)-(A5) of Theorem with T = oo. 

4.3. Proof of (Al). The operators Sf , T t N and G N are well defined on L 2 (S dN - 1 (y/£)) 
for any £ > thanks to the facts 



1 N f 

(G N ^<p) LH§dN ^ {VE)) = -- Jjr(K-^i) b( 



cos 9ij) [iplj -ip] 2 da 



and (Sf )* 



;<JV 



T t on this space, for any £ > (see for instance [Ml Chapter 9, 



't ) — D t 
section 2]). 

Then G N is well defined on the domain C\ 
using the fact that for any <p n G Dom(G N ), (p., 

that the restrictions of tp n , G N tp n to the sphere of energy £ belongs to L 2 on this compact 
space, and then, using the preceeding discussion, we deduce that ip = on this subspace. 
Since we can argue in this way for any £ > 0, we deduce that ip = 0, which is the definition 
of being closable. 



dN ), and it is closable on this space by 
->■ 0, G N (p n —> ip, one has for any £ > 
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Then it remains to prove bounds on the polynomial moments of the TV-particles system. 
We shall prove the following more general lemma: 

Lemma 4.2. Consider the collision kernel B = \v — b{9) with 7 = or 1 and b E 
^([0, 1), (1 - z) 2 ). This covers the three cases (HS), (tMM) and (GMM). 
Assume that the initial datum of the N -particle system satisfies: 

1 N 

supp/ " C [V e R Nd ; M 2 N (V) < £ } , M 2 N = ± £ \v 3 \ 2 

i=i 

and 

1 N 

</o Ar ,Mf><C 0ife <oo, = ^>2- 

i=i 

TTiere we /iawe 

sup (ff, Mjf ) < max{C ,fc; Sfc} 

where at G (0, 00) depends on k and £. 
Proof of Lemma \4-%[ Prom (|2.6p we have 

(4.7) supp /f C |V € M^; M 2 7V (y) < «?} , 

by taking := l^vf in (@3I). 

Next, we write the differential equality on k-th moment 

j t (/- jf E Kl fc ) = 7^2 E (/^ K - v h V K (v n ,v n) ) , 

with 

1 



K,{v h ,v j2 ) = - / 6(6» ilj2 ; 



1 * ifc 1 1 * \k 1 \k 1 1 

Kl + -\ v h\ -\ v h\ 



da. 



2 

From the so-called Povner's Lemma proved in [35\ Lemma 2.2] (valid for singular collision 
kernel as in our case), we have 

lC(v jlt v ja ) < Ci (lujil* -1 |« i2 | + |«i a |* _1 ) - Ci flu,-/ + \v j2 \ k ) 

for some constants C\, C2 £ (0, 00) depending only on k and b. 

By using the inequalities \vj 1 — vj 2 \ > \v 3l \ — \vj 2 \ and 1^ — Vj 2 \ > \vj 2 \ — {vjA in order 
to estimate the last term when 7 = 1, we then deduce 

\v h - v j2 \K,(v h ,v h ) < 
<Ca[{l + \u h \^ 1 )(l + \v ja \ a ) + {l + \v jl \ 2 ^ 

for a constant C3 depending on C\ and Ci. 
Using (symmetry hypothesis) that 

Vfc>0, (f t N ,\v 1 \ k ) = (f t N ,M»), 

and (|4.7j) we get 

I (f t N , \ Vl \ k ) < 2C7 3 </f , (1 + Mf +7 _,) (1 + M 2 ^)> - 2C7 2 (/f,M fe+7 ) 

< Ci (1 + £) (l + (jf , luil*^- 1 )) - 2 C 2 (if, M fe+ ^) . 
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Using finally Holder's inequality 



k-y+l 



<(/f> 



(fc-7+l)/(fc+ 7 ) 



we conclude that y(t) = (f i 



N 
t > 



\ v l\ 



satisfies a differential inequality of the following kind 



y' < -K x y dl + K 2 y e * + K 3 

with 9\ > 1 and 9 2 < 0\, which concludes the proof of the lemma. □ 

Lemma 14.21 proves (Al)-(i) with m e (v) = \v\ 2 , (Al)-(ii) with m\(v) = |u| 4 (and we do 
not need (Al)-(iii) in this particular case: we may take 711,3 = 0). 

4.4. Proof of (A2). Let us define 

PQ\ := {/ e Pi(JH d ) ; (/, m e ) < £ j endowed with the distance induced by | • I2. 
Let us prove (A2)-(i)-(ii) and (A2)-(iii) with: 
(4.8) V/ eP 6l , Vt 6(0,1] \S? L fo-f -tQ(f J )\ 2 <Ct 2 , 

(one can prove more generally that the application <S(/o) : [0, r) — >• Pg i; i h-» S(fo)(t) := 
S? L (f ) is C^C^rJjPoJ, with 5(/ )'(0) = Q(/o,/o) and r > 0). 

Let us recall the following result proved in |44j . We provide its proof for the the sake 
of completeness and because we will need to modify it in order to obtain similar results in 
the next sections. 



Lemma 4.3. For any /o,<7o 6 P2OR )> the associated solutions ft and gt to the Boltzmann 
equation for Maxwellian collision kernel satisfy 

(4.9) sup|/t- g t \ 2 < |/o-5o| 2 - 



Proof of Lemma \4-3\ We recall Bobylev's identity for Maxwellian collision kernel (cf [4]) 
T {Q + (f, 9)) (0 = Q + (P, G)(0 =: \ J b (a • tj [F + G- + F~ G+] da, 



with F = f, G = g, F ± = F(£±), G ± = G(£±), i = and 

r 



V + Kk), r = ^(e-iek)- 



2 v, „, „ , 2 
Denoting by D = g — f, S = g + f ', the following equation holds 



(4.10) 



d t D = Q(S,D) 



D+ S- D~ 5 H 
— = — + — — 



D 



da. 



We perform the following cutoff on the angular collision kernel: 



ok [o~ ■ £,) da = K, bx = b li 



\>5(K) 



for some well-chosen S(K), and we shall relax this assumption in the end (using uniqueness 
of measure solutions of [S]). Using that HSU 00 < 2, we deduce in distributional sense 



d \D\ 



\D\ 



\D\ 



^TT7Jo+K-^ < sup j-r=- 



dt\Z\ 



sup 



S d-1 



&RT ( (7 • £ 



+ 



do" 



with 



1/2 



1/2 
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d \D\ „ \D\ r , / \D\ 
~T7T9 + K T7T9 ^ K I SU P T7T2 



By using |£ + | 2 + |£ | 2 = 1, we deduce 

d \D\ 

from which we deduce 

for any value of the cutoff parameter K. Therefore by relaxing the cutoff K — > oo, we 
deduce □ 

Hence we deduce that SjF is C°' 1 (Pg 1 , P Ql ) and (A2)-(i) is proved. 

Lemma 4.4. For any f, g E Pi with same momentum and finite second moment, we have 

(4.H) |Q(/,/)li<C U(l + \v\)df(v)Y 

(4-12) \Q{fJ)v\ x <C ^ d (l + \v\)df(v^j 

(4-13) \Q(f, f)\ 2 <C (J (1 + \v\ 2 ) df(v)) 

and 
(4-14) 

\Q(f + 9,f-9)\i<C ( f (l + H) (df(v) + dg(v))) (\f-g\ 1+ [ (l + \v\)\df-dg\(v] 



(4.15) |Q(/ + <7,/-<?Hi<C (7 (i + M) (W + W)) (l/-5|i + |(/-5)Hi)- 

(4.16) |Q(/ + 5,/-5)l 2 <C (X d (1 + |u|) C^W+W)) (l/-5| 2 + |(/-5)Hi)- 



Proof of Lemma \4-3\ We prove the second inequalities (|4.14p - (|4.16p . The first inequalities 
(|4.1ip - (|4.13p are then a trivial consequence by using 

Qif, f) = Q(f, f) - Q(M, M) =Q(f-M,f + M) 

where M is the maxwellian distribution with same momentum and energy as /, and then 
applying P~T4|) or (^TU|) with f — M and f + M. 
We write in Fourier: 

T (Q(f + g,f-g)) = Q(D, S) = \ I b(a • |) {S(t) D(D + S(D ^(£ + ) - 2 £>(£)) 

where Q is the Fourier form the symmetrized collision operator Q, which yields 

Q(D,S) 

<Ti + T 2 + % 



with 



Tx < J ■ 6 \S(t + )\ ^^^rda<C\D\ 2 
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and 



r 2 < j b{a . a |5( y 21 Jii ^ < cp| a (£(i + ho + 



and 

Ts < 2 /- ^.a ^q^u 

■/§<*-! |?| 

which concludes the proof of (|4.14p . The proof of (|4.15p and (|4,16p are similar. □ 

The proof of (A2)-(ii) is a consequence of (|4.16j) . while the proof of (A2)-(iii) can be 
done by repeated use of Lemma l4,4i Indeed, we start from 



ft~k= f Q(fJ)ds 
Jo 



from which we deduce thanks to (|4.1ip , (|4.12p , (|4.13p 

\ft - /oli + \(ft ~ fo) u|i + \ft - fo\ 2 <Ct. 

Then write 

ft - fo -tQ(fo, f Q ) = f (Q{f s , f s ) - Q(/o, f Q )) ds = f Q(f s - f , f s + fo) ds, 

Jo Jo 

from which we deduce thanks to (I4.12p . (I4.13P 

\ft-fo-tQ(f Jo)\ 2 <C A|/ a -/o| 2 + |«(/»-/o)li) ds<C f sds<Ct 2 . 



4.5. Proof of (A3) with Q\. Define Ai(p) := (p, 1 + If] 4 ) 1 v for any p G Pg 1 and some 
i] G (0, 1). Let us prove that for any G (Pg ljr ;R), r G Rg 17 

(4.17) ||Mf (V^r 1 (G"**-7r"G~) *|| LOO(Ejv) < %^]^(P Sl , r;K )> 

for some constant C a and where Etv = {V G (R d ) N , \V\ 2 < £}. 
First, consider velocities v,v*,w,w* G M rf such that 

v + v* \v — vJ v + v* \v — vJ 

iu = a, w* = a, a G a 

2 2 ' 2 2 

Then 5 V + 5 Vt — 5 W — 5 Wt G ZPg x . Performing a Taylor expansion up to order two and one, 
we have 

e iv < _|_ e iv *'£ — e iw 't — e iw *'£ = 

= i(w- v)£ e m < + 0{\w - v\ 2 |e| 2 ) +»(«)♦- u*)^ e*** + 0(K - v*\ 2 |£| 2 ) 
= i(w-v) e iv < + 0(\w - v\ 2 \£\ 2 ) +i(w*- v*) (e iv -t + 0(\v - v.\ |£|) + 0(K - v*\ 2 |£| 2 ) 

= 0(\Z\ 2 \v-v*\ 2 cose) 

thanks to the impulsion conservation and the fact that 

\w — v\ + \w* — v*\ < V2\v — v* \ (1 — cos#). 
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We then deduce 

\5 V + S Vt -S w - 5 W J 2 = sup 

£<=R d 

Consider V G E^r and define 

,N 



< C\v — v*\ (1 — cos#). 



r V ■= ((/iy,Zl),. • • , (llv> z d), (A*V>k| 2 )) e Rf. 
Then for a given $ E C^ ,r? (P<j 1)IV ;R), we set := D$[/iy], Ujj = (vi — Vj) and we compute: 

1 ^ r 



M 4 



2JV 

=:7i(F) + / 2 (y). 
For the first term I\{V), thanks to Lemma 12.131 we have 

N 

h{V) = 2N^ S / + W) " ^ ~ Wi)] d ° 

i,j=i Jsd 1 

N i 



H'V* ~ H'V 



1+T) 



da 



2N 2 




6(0) [<t>{v*) + <t>(w*) - <f>(v) - 4>{w)} ^{dv)^(dw)da 



For the second term l2(V), using that 

M 4 (^ :j )<CM?(V), 

we deduce 



r N r r 

ij=l us 

< C[Mf(y)]^[^ (P6i , rv) ^(l + 0, 
which concludes the proof. 

4.6. Proof of (A4) in Fourier-based norms | • | 2 and | • [4. For /o,5o G P4(R d ), let us 
define the associated solutions /t and <^ to the nonlinear Boltzmann equation as well as 
ht := C^ L [fo] (go — fo) the solution to the linearized Boltzmann equation around ft- More 
precisely, we define 

dtft = Q(ft, ft), f\t=o = fo 
dt9t = Q{gt,9t), 9\t=o = 9o 

d t h t = Q(ft, h t ) + Q(ht, ft), h\ t=Q = h := g - f . 
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Lemma 4.5. There exists A G (0,oo) and for any rj G (0,1) there exists such that for 
an y fo-,90 G Pg 1 , r , r G Rg i; we have 

(4.18) |/t-0t| 2 < C v e~^ xt max(M 4 (/o),M4( 5o )) {1 ^ )/2 |/o-fldl2, 
and 

(4.19) \h t \ 2 < C v e~^ xt max(M 4 (/ ) 1 M 4 ( go ))( 1 -")/ 2 |/ - Soil- 

As a consequence, the operator C^ L defined by C^ L [fo](ho) := ht is such that C^ L [fo] G 
M 1 (Gi, r ,Gl,r)- 

Proof of Lemma \475[ Step 1. Estimate in \ ■ \±. We proceed in the spirit of [I2J [7]. With 
the notation M. = M4, A4 = M.&, introduced in Example \2.21\ let us define d := f — g, 
s := / + g and d := d — M.[d], and then D := -F(d), 5 := .F(s) and 

£> ■= F{d) = D - M[d}. 

The equation satisfies by D is 

(4.20) d tJ D = Q(D,S)-d t M[d] 

= Q(D,S) + {Q(M[d],S)-M[Q(d,s)}}. 

We infer from [32] that for any j G N rf , there exists some absolute coefficients (ay) 
depending on the collision kernel b such that 

/ b(co S 6) [(«<)' + («<)', - («*) - («*).] cfr = £ ay («»") («*-'")♦ • 

/<Hd— 1 * * 

We deduce that 

V \i\ < 3, V\M[Q(d, s)], =Q = M;[Q(d, s)] = a M M M M i-M 

j,\j\<\i\ 

together with 
V|i| <3, V\Q{M[d],S)^ =0 = Mi[Q(M[d],s)] 

= E a iJ M J [M[d]]M l ^[ s }= £ ayMj^M^M 
j.lil<|i| j>b'l<W 

since Mj[.M[d]] = M^[d] for any |i| < 3 by construction. As a consequence, we get 

(4.21) V£GR d , \MQ(d,s)]-Q(M[d],S)\<C\C\ A (^2\ M ^-9t]\)- 

\i\<3 

On the other hand, from [42\ Theorem 8.1] and its corollary, we know that there exists 
some constants C, A G (0, 00) such that 



(4.22) Vi>0 (J2 \Mi[ft-gt]\) <Ce' Ai (^ \Mi[f - g ] 

\i\<3 \i\<3 
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Gathering (|4.20p . (|4.2ip and (|4.22j) and performing the same cutoff on the angular collision 
kernel as in the proof of Lemma 14.31 we have 



dm)i +K m< sup 



dt w 



sup / b K ( a ■ £ 

Is*- 1 



+ 



da 



+ Ce- Xt \ M df°~9o] 
Let us compute (the supremum has been droped by the spherical invariance) 



■■= [ b K (a 



+ 



da = bx [ cr ■ £ 

! S d-i 



{l+(a-£ 



da, 



so that 



Xk-K 



[ b K ( 



da 



X 



S>00 J gd-1 

Thanks to Gronwall lemma, we get 



b [a-i 



(l-(a- 



da := -A G (-00, 0). 



sup 



!M)!W-*><[ S up!%P 

+ C 3 [j2\ M dfa-9o}\ 

\N<3 



-At 



K-\ K -\ K-X K -X 



Therefore, passing to the limit K — > 00 and choosing (without any restriction) A G (0, A), 
we obtain 



sup 



IA(0I 



< Ce 



-xt 



sup 



IA>(0I 



+ J2\M t [fo-9o}\ 

\i\<3 



l£l 4 \ sm d l£l 4 

or equivalently (and with the notations of Example 12.2 lj) , 

||KH|4<C7e- A '|||do|||4. 
Step 2. From the preceding step and a trivial interpolation argument, we have 

\f-g\ 2 <\f-9-M[f-g]\ 2 + C l^2\Mi\f-g]\ 

Vni<3 



<Wf-9-M\f-gW£ \f-g-M[f-g]\T + C ^|M,[/- 5 ]| 

\l*l<3 

< C (1 + M 4 (/ ) + M 4 ( 50 )) e" (A/2) *. 
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We conclude by writing 

\f-g\ 2 <\f-g\l\f-g\\-\ 

using Lemma 14.91 for the first term and the previous decay estimates for the second term. 

Step 3. The same computations imply at least formally the same estimate on ht as stated 
in Lemma 14.51 Now, we proceed to the rigorous statement. We consider a truncated 
(cut-off) model associated to b £ = bl cos 0<i_ £ G L (5 1 1 ) and an initial datum G 
L3,(R 3 ) such that (M, 1) = 0. By standard argument there exists a unique solution h\ G 
C([0, oo); L^M 3 )) to the associated linearized Boltzmann equation arround ft, and this 
one fulfills (I4.19|) . Letting e goes to 0, we get that for some ht G C([0, oo); ^(M 3 ) — w) there 
holds h\ — ht in 5'(R 3 ) and fat is a weak solution to the linearized Boltzmann equation 
(corresponding now to b) and which satisfies again (I4.19D . In other words, we have built 
a linear flow Cf [fo] ■ Qi — > Gi si that for any ho G Q\ the distribution C^ L [fo](ho) ■= h t 
is the solution to linearized Boltzmann equation. □ 

Lemma 4.6. There exists A G (0, oo) and for any r] G (0, 1) there exists C v such that for 
any f ,g G Pg 1)T) r G Rg ± , we have 

\ut\4 < O e | go - 7o| 2 

:= 9t- ft~h t = S? L (g ) - Sf L (/ ) - £f L [/o](5o - /o) 
('as proved below, 0Jt always has vanishing moments up to order 3). 

Proof of Lemma \4-6\ The following arguments can be fully justified as in |44j by truncating 
b and passing to the limit. Consider the error term 

to := g — f — h, £1 := to. 

Again we perform the angular cutoff of Lemma 14.31 with cutoff parameter K, then the 
evolution equation (in the Fourier side) satisfied by f2 is 

(4.23) d t n = Q(n,S) + Q + (H,D). 

Let us prove that 

< 3, Vt > 0, Mi[oj t ] ■= / v { duj t (v) = 0. 

J~BL d 

Consider the equation on to t : 

d t u t = Q(u J ,(f + g)) + Q(h,(f-g)) 

and the fact that, for maxwell molecules, the i-th moment of Q(fi, fz) is a sum of terms 
given by product of moments of f\ and fi whose orders sum to Hence using that for 
some absolute coefficients (aj j) we have 

i b(co S 0) [(v i )' + (v% - (v l ) - (v%] da=J2 <kj H (w <_i ), , 

i>\i\<i 



we deduce 

3 

dt 



V \i\ < 3, -^M< [wt] = ^ oij Mj [to t ] M^j [f t + g t ] + ^ Mf [/i t ] Mj_j [/ t - 

3<i j<i 



and since 

V|i|<l, Mi[ht] =M i [f t -gt] = 
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we deduce 



V \i\ < 3, ^Afi[w t ] = a iJ M M] Mi-jtft + 



9t\ 



3<i 



and from the initial data ujq = this concludes the proof. 
We now consider the equation in Fourier form 

d t n = Q(n,s) + Q + (H,D) 

and we deduce in distributional sense 



d\n(0\ K mt)\\ <Ti , T 
Jt ^r + K^ r )<T 1 + T 2 , 



where 



7i := sup 



b [*.£ 



3 Jsd-1 



0(^+)5(t 



+ 



da 



sup 



< 



< sup 



b [a-i 



+ 



ie+r ler ie-r lei 



TOP 

ici 2 



sup 



6 (o--^ 



+ 



da 



da 



where A^- was defined in Lemma 14731 and 



fT(r)r>(r)+^(r)£>(e + ) ^ 



i 

< o sup 



< 1/itU Mil 



6 



mm in 2 to+)i 2 |g(ni 2 in 2 ', (/rT 



|£+|2 |£-|2 |£|2 |£+|2 |£-|2 |£|2 

& ( f • Co ) ( 1 - O" • |o) ^(7 



Cg-MAt |^ nU | dn | 3 < c-e-a-^At !do |i+fl 



l 0|2 I "0| 2 

by using the estimates of Lemma 14.31 
Hence we obtain 



* lei 



lei 4 ) 



lei 4 ; 



We then from the Gronwall inequality and relaxing the cutoff parameter K as in Lemma 
31 (assuming without restriction (1 — rj)\ < A 



S e |5q — Job 



□ 



if) 
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4.7. Proof of (A5) in Wasserstein distance. We know from [32] that 

S upW 2 (S? L fo,SF L go) < W 2 (f ,g ). 

As a consequence, from Wi < W 2 < C k wf k) and Wi = [■]*, we deduce that (A5) holds 
with Q(x) = x^ k \ Tz = Lip and Pg 3 = P 2 endowed with || • ||g 3 = [■]*. 

4.8. Proof of (A5) in negative Sobolev norms. It is also possible to prove easily, in 
the cutoff case, that the weak stability holds in Sobolev space on finite time: 

Lemma 4.7. For any T > and ft, gt solutions of the Boltzmann equation with maxwellian 
kernel and initial data fo and go, that there exists s G (d/2,d/2 + 1) and a constant Ck 
such that 

sup - 9t\\ H — ^ C T,s ll/o - 9o\\h— ■ 

[0,T] 

Remark that this proves (A5) with Q3 = P(M. d ), endowed with the H~ s norm. This is 
useful in order to obtain the optimal rate 1/y/N of the law of large numbers. 
Let us only sketch the proof. We integrate (|4.10p against D / (1 + \£\ 2 ) s : 

d„ ,„ Iff ( A \D-S + D + D+S- D-2\D\ 2 ] 
and we use Young's inequality together with the bounds 

IML. II^1L<II/+5||m 1 <2 

to conclude. 

5. Hard spheres 

5.1. The model. The limit equation was introduced in Subsection 11.11 and the stochastic 
model has discussed Subsection 14.11 We consider here the case of the Master equation 
(14. 2p . (|4.3p and the limit nonlinear homogeneous Boltzmann equation (jl.2p . (jl.3p . (|1.4p 
with T(z) = \z\. 

5.2. Statement of the result. We have the following theorem: 

Theorem 5.1 (Hard spheres homogeneous Boltzmann equation). Let us consider /o £ 
P(R d ) with compact support. We denote by £ and A two positive constants such that 

M 2 (f )<£ and supp f C {v G R d , \v\ < A}. 

Let us also consider the hierarchy of N -particle distributions = S^f^) issued from 
the tensorial initial datum f$ = f$ N , N > 1. Let us finally fix some 5 G (0, 1). 

Then there are some constants k\ > 2, C > 0, only depending on 5 and £, some constant 
C v depending on rj, and some constants a > depending on the collision kernel such that 
for any 

^ G N*, ^ = W 8p2®'"®WG^ 1 ' 0O (lT. 

we have 

(5.1) VN>2£, sup ((sFVft-iSFVo)) 9 "),? 

[0,T] v v ' 



< ce\\(p\ 



||/o]Imi II/oIIm, 1 . 
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and the dominant error is the last term given by 



e 3 (N) = e 

for any rj > small (the constant blows up when r] —> 0). As a consequence we deduce 
propagation of chaos with explicit rate. 

In order to prove Theorem l5.14 we shall prove each assumption (A1)-(A2)-(A3)-(A4)- 
(A5) of Theorem 13,271 step by step. Its application then exactly gives Theorem 15.11 We 
fix F 1= F 2 = C b (R d ) and Jjj = Lip(M d ). 

5.3. Proof of (Al). In the proof of Lemma 14.21 we have already proved that 

Vt > 0, supp/f C E N := <|V G R Nd ; M^(V) < «?} , 

which is precisely (Al)-(i) with 7772(7;) := \ v \ 2 -, as we h as that for any k > 2, 

su P (/f,Mf)<Cf 

i>0 

where depends on k, £, on the collision kernel and on the initial value {f^M^) 
which is uniformly bounded in TV (in function of k and A for instance). That is precisely 
(Al)-(ii) with 777.1(7;) := |7j| fcl where k\ > 2 will be chosen (large enough) in section [5771 
As for (Al)-(iii), we remark that for a given N-particle velocity V = (V\, V/v) G M. dN , 
we have 

Vesuppf® N <=> Vi = l,...,N, y.Gsupp/o Vi = l,...,yV, m 3 (y i )<m 3 (i), 
with 777,3(7;) := e ' 11 '. We conclude that 

supp/ ®" C |v G 1^; M^(y) < 7773(A)} , 

and (Al)-(iii) holds. 

5.4. Proof of (A2). We define 

P Sl := {f e P kl (R d ); M 2 (f)<8) 
that we endow with the total variation norm || • \\g 1 := || • ||^i. 

• The proof of assertion (i), that is for any a > a kl > and for any t > 0, the application 
/o 1 y S^ L fo maps BPg lta := {/ G Pg x , M^if) < a} continuously into itself, is postponed 
to section 15.61 where it is proved in (I5.4p a Holder continuity of the flow. 

• For any /, g G Pg x we have 

\\Q(g,g)-Q(fJ)\\AP = \\Q(g - f,g + f)\\ M ^ 

< 2 / / / b(8) I v — v*\ |/ — g I I f*+ g*\ da dv* dv 

JR d JK d JS^ 1 

(5-2) < 2(l + £)\\b\\ L i\\(f-g)(v)\\ Ml . 

We deduce that 

\\Q(g, g) - Q(f, f)Un < 2 (l + £?' 2 \\b\\ L , \\f - g f£ 

which yields Q G C^iPg^Qx). 
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Finally, we have for any /o £ Pg 1 



Wft - IqWmI 
and then, using (|5.2p . 

\\ft — fo — tQ(fo, fo)\\ M i 



[ t Q(fsJs)ds 


<- / 


Jo 


M{ J0 



< f \\QUsJs)\\ M i ds < 2,t\\b\\ L ,{l + Sf 

!} JO 1 



(Q(fs,fs)-Q(fo,fo))ds 

< [ t \\Q(fsJs)-Q(fo,fo)\\ M ids 
Jo 

< [ 2(l+£)\\b\\ L i \\f s -f \\ M xds<3t 2 
Jo 



\h (l+£) s 



which precisely means that t ' — y Sf L (/ ) is C 1 ' 1 at t = 0+ in the space Pg 1 . A similar 
argument would yield the C ' ([0, T],Pg t ) regularity on a small time interval [0,r]. 

5.5. Proof of (A3). For any /c G N*, any energy bound £ (which has been fixed once for 
all in the statement of the main result), any energy r d+ \ € [0, £] and any mean velocity 
(n,...,r d ) e B Rd (0,r d+1 ) we define 

POuv ■= {/ G Pk(^ d ); (f,v j )=r j ,j = l,...,d, M 2 (f) = r d+1 }, r := (n, r d+1 ). 

We also denote by the set of all admissible vectors r 6 M. d+1 constructed as above. We 
claim that there exists C = C^g such that for any r/ £ (0, 1) and any function 



with A(/) := M kl (f), we have 

(5.3) VVeEiv, |G 7V ($o M ^)-(G 0O $)(^)| < C I sup [$] r i,„ 



K<v) 



which is precisely assumption (A3) with e(N) = C/iV. 
Consider V G and define 

r v := ((A*v, Zl), • • • , Zrf), (//y, |^| 2 )) e R £ . 

Then for a given $ £ C^ ,r? (Pg lirv ; M), we set 4> := D$[fj,y] and we compute: 

i N r 



+ 



^ E 1^ - ^-1 - <t>) b(fiij) da 

jv . 
V \Vi-Vj\ / max \M kl 
',.7=1 Jb 



2N 

Wd?^) N 



M, 



a {$)} 



O 



M 1 



--:h(V) + I 2 (V). 
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For the first term Ii(V), thanks to Lemma 12.131 we have 

N 

h{V) = Wt ? lVt ~ Vjl J d -i h( - 9ij) [<PW) + HV*)-HVi)-m)} da 

' Nfj..\ ,.N 



2N 2 



d-l 



v-w\ I b{9) [<f>(v*) + (f>(w*) - <f>(y) - <t>(w)} $ (dv) (dw) da 



For the second term l2(V), using that 



M kl (^) := < (V£) = 1 \V n \ k ^j + \V*\^ + \V*\ k ^j 

nl+fel/2 / \ 



we deduce 



/ 1 JL / 4 \ 1+ A 

|/ 2 (F)| < cr<(v)[*] oi , (w ^Ej^-^l(]vJ j 



i i * 



< cr<(vo [*] ,, (JW ^ ^ E TO + 



< cM^vm cl , niPgirv) ±(i+£). 

We conclude that (|5.3p holds by gathering these two estimate. 



5.6. Proof of (A4) on a finite time interval [0, T], 

Lemma 5.2. For any given energy £ > and any 5 > i/iere exists some constants ki > 2 
(depending on £ and 5) and C (depending on £), such that for any /o,<?o £ -P<7i0^ d )> we 



(5.4) ||5f L (5o)-^(/o)|| Ml <e c W y/max{M kl (f ),M kl (g )} \\fo - 9o\\ 

(5.5) ||£Sf L [/o] (fo-g )\\ M1 <e c ^ sj MM) ||/o - sofc*, 



1-5 
Af 1 ' 



(5.6) 115^(50) - 5f L (/ ) - ^[/oKffo - /o)| 



M 1 



,2-8 

Im 1 



< e C{1+t) ^m a x{M kl (f ),M kl (g )} ||/ - go||] 
Proof of Lemma 15.21 We proceed in several steps. Let us define 

V/eM 1 ^), \\f\\ M y-= I (v) k \f\(dv), \\f\\ Ml := [ (v) k (l + log(v)Y\f\(dv). 
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Step 1 . The strategy. Let us define 

dtft = Q{ft, ft), f\t=o = fo 

9t9t = Q(9t,9t), 9\t=o = 9o 

dth = Q(ft, h t ) + Q(h t , ft), h\ t=0 = h := g - f . 

Existence and uniqueness for ft, gt and ht is a consequence of the following important 
stability argument that we use several times. This estimate is due to DiBlasio [13] in a 
L 1 framework, and it has been recently extended to a measure framework in [141 Lemma 
3.2]. Let us recall the argument for h. We first write 

(5.7) | j(v) 2 \h t \(dv) < J!] \h t \(dv)f t (dv*)\u\b(e)[(v') 2 +(vi) 2 -(v) 2 -(v*) 2 

+2 \h t \{dv) f t {dv*)\u\b{0) (v*) 2 dadvdv* 



da 



(this formal computation can be justified by a regularization proceedure, we refer to |14| 
for instance). Because the first term vanishes, we conclude with 

(5-8) j t \\ht\\ Aq <C\\f\\ Aq \\h t \\ M i. 

When ||/t||_kfi G L 1 (0,T), we may integrate this differential inequality and we deduce that 
h is unique. 

More precisely, we have established 

(5-9) sup \\h t \\ M i < \\g - f \\ M i ex.pl C / ||/ s || M i ds ) , 

[o,T] 2 2 V Jo 3 / 

and similar arguments imply 

(5.10) sup \\ft - g t \\ M i < \\go - fo\\ M i exp ( C / ||/ s + g s \\ M i ds ) . 

[o,T] 2 2 \ Jo 3 / 

It is worth mentioning that one cannot prove H/tH^i G tL 1 (0, T) under the sole assumption 
II/oIIm 1 < 00 because since it would contradict the non-uniqueness result of [30]. However, 
thanks to Povzner's inequality, one may show that ||/t||jvfi G L 1 (0,T) whenever H/oIIm^j 
is finite, with the definition 



II/oIIm^ := J (v) k log ((«))' df (v) < +oo 

(see ()5.13p below or [351 EZ]), which will be the key step for establishing (I5.4p and (15.5 
Now, our goal is to estimate (in terms of \\go — /oIIm 1 ) the M 1 norm of 

Ct ■= ft- gt- ht. 
The measure Ct satisfies the following evolution equation: 

dtCt = Q(ft,ft) - Q(gt,gt) - Q(h,ft) - Q(ft,h), Co = 0. 
We can rewrite this equation as 

dt(t = Qttt,ft + gt) + Q(ht,ft-gt). 

The same arguments as in (|5,7p - (|5.8p yield the following differential inequality 

j t ||Ctll M i < c KtWui Wft + gthq + QihJ't - gt) Ml , IICoIIm, 1 = °- 
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We deduce 



(5.11) 


sup WCtWm 




te[o,T] 


Since 












(5.12) 


f 


Q(h s ,f s ~ 




Jo 





Q(h s , f s - g s ) ^ ds J exp ( C 



+ 9s\\mI ds 



Ml 



ds < C I sup ||/i t || M i 
\[0,T] 2 J 

+ C sup \\g t - /t|| M i 

\[0,T] 2 



Wfs ~ 9s\\mI ds 



ds 



WfsWui + II^IImi 



we deduce from (|5.9|) and (|5.1U|) 

SU P IICt|| M i ^ c llfo - Mm] ex P ( [ 

[0,T] 2 2 \J0 

{ {[ II* ~ /.II**) exp ||/.||^) + [[ exp 

Hence the problem is reduced to time integral controls over H/tH^i, H^Hm 1 ? Wft ~ 9t Hm 1 
and ll/ttll^i. 



II/sIImI + \\9s\\ M £ 



til M3 1 

Step 2. Time integral control of f and g in Mg . In this step we prove 
(5.13) 



J \\ft\\ M l^dt<C £ T + C'\\h\\ Mii £=1,2, 

for /, where Cs is some energy dependent constant, and C is a numerical constant. The 
same estimate obviously holds for g. These estimates are a consequence of the accurate 
version of the Povzner inequality as one can find in |35[ [27] . Indeed it has been proved in 
[351 Lemma 2.2] that for any convex function \£ : R rf — > R, ty(v) = ip(\v\ 2 ), the solution f t 
to the hard spheres Boltzmann equation satisfies 

' ' " 11 f t (dv) f t (dv*)\v -v*\Kv(v,v*) 



dt 



with K^j = G^j — H^,, where the term "behaves mildly" and the term Hq, is given by 
(see f35 :) formula (2.7)]) 

ctt/2 



ip{\v\ 2 cos 2 9 + \v*\ 2 sm 2 6) -cos^V(M ) -sin 2 #V(H ) 



dO, 



(note that Hq, > since its integrand is nonnegative from the convexity of ip). More 
precisely in the cases that we are interested with, namely if?(v) = ip2,i(\v\ 2 ) with ipk/ir) = 
r k l 2 (log rY and £ = 1,2, it has been established that (with obvious notations) 



(5.14) 



|<W«,«.)| < Ci ( v ) (Ml + ( V ) 2 )Y M (Ml + M 2 )) 1 



On the other hand, in the case £ = 1 we easily compute (with the notation x := cos 2 8 
and u = \vA/\v\) 



Vx,«Gl, x G [1/4,3/4], u G [0,1/2], 



^2,1 ( | f | 2 cos 2 + |t>*| 2 sin 2 • 



cos 



0*h,i {\v\ 2 ) - sin 2 0^ 2 ,i (H 2 ) 



(1 - X) ^2,1 {u 2 ) + X ^2,1 (1) - lh,l {{1-X)U 2 + X) 



> Kq \v\ 
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for some numerical constant hq > (note that this lower bound only depends on the 
strict convexity of the real function ^2,1 )• We straightforwardly deduce that there exists 
a numerical constant k\ > such that 

H 2 ,l(v,V*) > K\ \v\ 2 l|„|>2|t,„|. 

Similarly, in the case I = 2, we have 
Vx,aeK, x G [1/4,3/4], u G [0,1/2], 

^2,2 (|^| 2 cos 2 + |^ + | 2 sin 2 6>) -cos 2 6*^2,2 (M 2 ) - sin 2 9 ^2,2 (|u*| 2 ) = 



+ M 2 



2|u| 2 log |u| 2 {(1 - X)lp2,l {U 2 ) +X^2,l(l) -^2,1 ((1 ~X)U 2 + x)} 

(1 - x) V>2,2 (u 2 ) + x ^2,2(1) - V"2,2 ((1 - x) n 2 + x) 



> 2 kq 1 17 1 2 log \v\ 2 , 



and then 

H2,2(V,V*) > 2/Cj \v\ 2 log|v| 2 l|„|>2|„,|. 

Putting together the estimates obtained on G 2 ,t and we deduce the Povzner in- 
equality 

(5.15) \v - v*\K 2/ < C {{v) 2 + (v*) 2 ) - k \v\ 3 (log(v)Y~\ 
and we finally obtain the differential inequality 

j t \\ft\\ Mll <2C(l + £)-KM 3/ _ 1 , 

from which (|5.13p follows. 

Step 3. Exponential time integral control of f and g in M3 (proof of (|5.4j) and (|5.5p ). Let 
us first prove that 

(5.16) e° / ° T|i/sll «3 1 ds < ^MUfo), e° ^ M "i ds < y^gV) 

for any k > kg, with kg big enough. 

This is a straightforward consequence of the previous step and the following interpola- 
tion argument. For any given probability measure / G Pk(M. d ) with M 2 (f) < £, we have 
for any a > 2 



M i = / <w> 2 (l + log«7;> 2 )) (l ( , )2 < a + l ( „ >2 > a ) f(dv) 

< (l + f)(l + logo) + - / (v) 4 (l + log((v) 2 )) f(dv) 

< (l + £)(l + loga) + i||/|| Ml . 

a 

By choosing a := ||p|| M i, we get 

(5-17) IHU^ <2(l + f) (l + log||p|| M ^ 

On the other hand, the following elementary Holder inequality holds 

(5.18) Vk,h> G N, k' < k, V/ G MjJ, ||/|| M i < U/fc* 7 * ||p||$. 

Then estimate (foTTBT) follows from fjoTTH]) . (l5T7jl and (jSTTgl) with fe 7 = 2 and & = k £ > 5 
large enough in such a way that 

c'2{i + e)\<\, 
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where C is the constant which appears in (I5.13p . We then deduce (15. 4p from (|5. lOj) . and 
(similarly) (|5.5h from (|5.9p . 

Step 4. Time integral control on d and h. Let us prove 



(5.19) J (\\dthq + W h thq) dt 



. C e (l+T) C '[y^ul+\\9o\\ M i ) / \ 

for some energy dependent constant and some numerical constant C . Performing 
similar computations to those leading to (|5.7|) . we obtain 

d 



^iNIm*! ^ / j \h\(dv) f t (dv*) \u\K 2 ,i(v,v*) 

+2C Iff \h t \(dv)f t (dv*)\u\{v*) 2 (l + log^.) 1 



Thanks to the Povzner inequality (j5. 15j) . we deduce for some numerical constants C, n > 
Integrating that differential inequality yields 



IMIa/i, + « / INU/i * < C sup \\h t \\ M i 

\[0,T] 2 



Using the previous pointwise control on H^Hmi an d Q5.13P (with £ = 2) we deduce 
(5-20) T tft < ( 1+T ) e Cl " /o Hi ||/ || M i 2 ||<? - /o|| M i + ||<7o - /oll^ • 



Arguing similarly for dt, we deduce f|5. 19|) . 

Step 5. Conclusion. By gathering the estimates (|5.1ip - (j5.12p - (j5.16]) - (|5.19p . we obtain 

II/- II ^ C £ (l+T) °' [\\Mm1 + \\9o\\ m i 
su PllCt|lA/i < e° £l + ; e V 2,1 2,1 
[0,21 

(||/oIIm! i2 + [Iso[Imi >2 ) ll^o - /oIIm| ho - MIm^ > 

for some energy dependent constant Cg and some numerical constant C . Then arguing 
as in the end of step 3, using (|5.17p and f|5. 18j) with k large enough, we get 

sup ||C t || M i < e c ^ Vmax{M fc (/ ),M fc ( 5o )} ||/ - ^oll^^ » 
[0,21 

from which estimate (15.61) follows. □ 



5.7. Proof of (A4) uniformly in time. Let us start from an auxiliary result. It was 
proved in |36| that the nonlinear and linearized Boltzmann semigroups for hard spheres 
satisfy 

(5.21) llsf^Li, -i, < C z e~ xt , lle^^L^ < e~ Ai 

v ' II t \\L L (m z j — z 'II \\L L (m z j — z 

where m z (v) := e 2 ' 1 '', z > 0, A = X(£) is the optimal rate, given by the first non-zero 
eigenvalue of the linearized operator C in the smaller space L 2 (M _1 ) where M is the 
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maxwellian equilibrium (see [36, Theorem 1.2]), and C z is some constant depending on z 
and the energy £ . 

Lemma 5.3. For any given energy £ > 0, there exists some constants k\>2 (depending 
on £ and 5) and C (depending on £) and n 6 (0,1), such that for any fo,go £ Pg 1 (M. d ) 
satisfying 

Vi = l,...,d, (f ,Vi) = (g ,Vi) and </ , \v | 2 } = (g , \v\ 2 ) < £, 

we have 

(5.22) \\S? L (g ) - S? L (f )\\ Ml < e C ~^ Jmax {M fcl (/ ), M kl (g )} \\g - f \\£ , 



1+2 

(5.23) \\£Sr[M(9o ~ /o)|| M i < e°-t* ^M fcl (/ ) \\g - /oll A/ 2 i , 

(5.24) \\S t NL (g ) - S t NL (f ) - CS^ L [f ](g - /o)|| Ml 

< e c -^y m a X {M fcl (/ ), M kl (g )} \\g - / ||jJ?. 
Note that ([5T241 implies (A4) with T = oo, Pg 2 = Pg x . 

From [I], there exists some constants z, Z (only depending on the collision kernel) such 
that 



Bvp\\ft + gt + ht\\ L i <Z, m z (v):=e 
t>i 2z 



We also know from (|5.2ip that (possibly increasing Z) 

V* > 1 II/* " M||^ +\\g t - M\\ Lln2z <2Ze-", 
where M := Mj = M go stands for the normalized Maxwellian associated to /o and go and 
(5.25) H e£ *|L <Ce~ xt , Ch:=2Q(h,M). 

We write 

dt(ft ~ 9t) = Q(ft ~ 9t, ft + 9t) = C(ft ~ 9t) + QUt ~ 9t, ft~M) + Q(f t - g t ,g t - M) 
and we deduce for 

u(t) := \\ft-9t\\ M }, 
the following differential inequality (starting at some time To) 

u(t) < e- X ^- T °K(T )+C fe- A ^ \\Q(f s - g„ f s - M) + Q(f s - g s ,g s - M)\\ M1{mz) ds 

J To 

(this formal inequality and next ones can easily be justified rigorously by a regularizing 
proceedure and using a uniqueness result for measure solutions such as [T71 [2] ) . Therefore 
we obtain 

u(t) < e -Mt-Zb) U (T )+C f e" A M ||(/ s - M) + (g s - M)\\ Ml{{v)mx) \\f s - g s \\ M i ({v) mz) ds. 

By using the control of M l ((v) m z ) by M 1 (m2 Z ), the decay control (|5.21|) and the trivial 
estimate 

-A(t-To) < -A( t _T ) -A( s _t ) 



we get 



«(t) < e-M*-To) u(To) +Ce -|(t-To) r e -|( s -T ) ||/. 

•/T 
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We then use the following control for any a > 0: 

\f-g\ (v)e*W 



9\\m} s 

(v) m z 



<a[ \f-g\ e ^\ + e- za [ (/ + g) 

J\v\<a J\v\>a 

< au + e~ za Z. 



Hence we get 



C u + e z Z < (1 + Z) u when u > 1, (choosing a := 1) 

\\f-g\\ M l < { r) , , , , 

wm* ^ i i i Q g u \ u -\- u z when u < 1 (choosing — |a := log u) 
and we deduce 

||/-<?|| M i <Ku (l + (logu)_), tf:=l + - + Z. 

Then for any 5 > small, we conclude, by choosing To large enough, with the following 
integral inequality 



(5.26) u(t) <e- A ( t - To )u(r ) + 5e-^*- To ) / e-t( s -!) ns (l + (logu s )_) (is. 

JTq 

Let us prove that this integral inequality implies 

(5.27) Vt>l, u(t) < Ce-i'ufTo) 1 - 5 . 

Consider the case of equality in (|5.26p . Then we have u(t) > e~ A C -1 ) u(l) and therefore 
(l + (logu a )_) < (l + (Iog«(r ))_ + A(s-T )). 

We then have 

tf(t) := T e~* ( s "^o) Us h + (log Us )_) ds < f e"* u s (l + (log«(r ))_ + A (s - T 

< (l + 0ogu(T o ))_) / e-^^-^Usds 
Jto 

and we conlude the proof of the claimed inequality (|5.27p by a Gronwall-like argument. 

Then estimate (|5.22p follows by choosing 5 small enough (in relation to rj) and then 
connecting the last estimate (|5.27p from time To on together with the previous finite time 
estimate (|5.4p from time until time To. 

Then estimates (|5.23p and (|5.24|) are proved in the same way by using the equations 

d t h t = C(f t -g t ) + Q(h t ,f t -M) 
(which is even simpler than the equation for ft — gt) and 

OtCt = 2C(h t ) + Q((t,ft - M) + Q((t,g t - M) + Q(h t ,d t ). 
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5.8. Proof of (A5) uniformly in time. Let us prove that for any z, Ms £ (0, oo) 

there exists some continuous function G : M + — >• R + , ©(0) = 0, such that for any fo, 
go £ PmzO^- d ), with rriz defined in the previous section, such that 

rn z IIL z 

there holds 

(5.28) supWj < 9(Wo), 

t>o 

where W% stands for the Kantorovich-Rubinstein distance 

W t = W 1 {S? L (f ),Sj > L (g )). 
As we shall see, we may choose 

(5.29) Q(w) := Gminjl, e -( lo g( w / e o))_ /2 J , 

for some constants 0, 6q > (only depending on z and Aig). 
We start with 

(5.30) Vt>0 W t < \\(ft-9t)\v\\\ M i < l\\(ft + 9t)(v) 2 \\ M1 = 1 + 8 =:Q. 

Let us improve this inequality for small value of Wo. On the one hand, it has been 
proved in |17l Theorem 2.2 and Corollary 2.3] that 

(5.31) W t < W + K fw s (l + (logW 8 )_) ds, 

Jo 

for some constant K. Whenever Wt < 1/2, the integral inequality fl5.31f) implies (possibly 
increasing the constant K) 

W t <W + K f W s (log W s )_ ds. 
Jo 

From the Gronwall lemma we deduce 

(5.32) W t < (W ) cxp( ' Kt) whenever W t < 1/2. 

On the other hand, from [5] and [36], there exists A2, Z > 0, z G (0, z) such that 

Vt > ||/ t - M h \\ L}nz + \\g t - M go \\ LLz < Ze- X \ 

where Mf and M go stand again for the normalized Maxwellian associated to /o and go. 
Denoting by Uf and u go the mean velocity of /o and go, by 9f and 9 go the temperature 
associated to /o and go, and by £f and £ go the energy associated to fo and go, there also 
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holds 



W 1 (M f0 ,M l 



.go ) 



< C s (\\M fo -M go f H . s y 
id 2 



< C, 



< C, 



g "90 2 "90 V g 



' U 90 V"90 5 



(0 



2s 



2 \ »7 



/ 



<?/o-^oriei 4 + i« /o v^b 



'|2 |t|2 



*3o V 9o I 



(0 



2s 



< C s {\9 fo -9 go \*> + \u fo 

< C s {\£ h -£ 9 S + \u h 



u 



"go I 



127) 
1 27) 



) 



) 

< Cs (W 2 (/o,3o) 2r? + W r i(/o,5o) 2 ' 7 ) 

< C s (w 1 (/ , 9o ) ,?/2 + W 1 (/o,3o) 2,? 



(see also [12] for more general estimates of the Wasserstein distance between two gaus- 
sians). Gathering these two estimates, we deduce 

(5.33) Vt > W t < Z 1 e~ xt + Z 2 W% /2 . 

Let us (implicitly) define T, W E (0, ^) by 

Z ie ^ AT = i and next (W ) exp( ^ T) = ^. 
Then, for any Wo € (0, Wo) we have 

Let us search for t* such that the two functions involved in the minus function are equal: 
<j>(t*) := | log Wb | exp(-Kt*) = At* -logZ=: ^(t*). 



Vt > W t < minlo^r^ ; ^ e^'j + Z 2 W^ 2 



The time i* is unique because 4> is decreasing while V is increasing (and it exists for 
log Wo < logZ). Choosing 

t* :=llog(|logW | 1/2 ), 

we find 

0(t«) = |logW | 1/2 > A log (|logW | 1/2 ) -logZ = V(t tt ), 

at least whenever Wo E (0, Wq], with Wq E (0, Wo] small enough. As a consequence, for 
any Wo E (0, Wq] we have 

Vt>0, W 4 < 20(t tt ) = 2e- |logH/o|1/2 . 
This concludes the proof. 
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6. Extensions and complements 

In this section first we generalize the chaos propagation results on the Boltzmann 
(maxwellian or hard spheres molecules) to the case where the limit 1-particle distribu- 
tion is not compactly supported. This shall rely on a construction due to Kac of the 
iV-particle initial data together with careful study of the dependency of the constants in 
terms of moments of the data. Second, as a corollary, we use our global in time results to 
give a new method for studying chaotic convergence of the iV-particle equilibria towards 
the limit 1-particle equilibrium. The old question of connecting the long-time behavior was 
raised by Kac and it motivated its whole study of chaos propagation for particle systems. 
In the case of classical gas dynamics, we thus recover a well-known computational results 
(namely the marginals of the constant probabilities on y/~NS Nd ~ 1 converges to products of 
gaussians) without any explicit computations, only using the properties of the Boltzmann 
semigroups. This new method shall prove highly useful in the context of granular gases 
where the steady states or homogeneous cooling states are not explicitly known. 

6.1. Construction of chaotic initial data /q with prescribed energy. For the sake 
of completeness, let us recall, following [22J, how to construct a /o- cna °tic sequence of 
initial data /q (i.e., which has the "Boltzmann's property" in the words of Kac). 

Lemma 6.4. Consider fo G P(IR d ) with finite energy M me (/o) := M2(/o) = £ G (0, oo) 
and which fulfills the following moments conditions M mi (/o) = M^^. < oo, % = 0,1,3, 
for some radially symmetric and increasing weight functions mi, i = 1,3, and mo(x) := 
exp(a|x| 2 ), a > 0. Then for any given increasing sequence (aj\r)N>i (which increases as 
slow as we wish in general and may be chosen constant when fo has compact support), 
there exists a sequence /q G P(W dN ), N > 1, such that 

(i) The sequence (f^)N>i is fo-chaotic. 

(ii) Its support satisfies 

supp/<f C S^- 1 (yNS) := {y G R dN ; M?(V) = s} C E N . 
(hi) It satisfies the following integral moment bound based on mi : 

ViV G N*, (f N ,M^)<Cst{M N ^). 
(iv) It satisfies the following moment bound on the support involving m^: 

supp/o^ C K N := [V G R dN ; M% 3 (V) < a N } . 



Sketch of the proof of Lemma \6.4\ We essentially recall briefly the key arguments presented 
in [22, Section 5 "Distributions having Boltzmann's property"] and check that the moment 
conditions required in the sequel of our paper can be satisfied. For the sake of simplicity, 
we assume with no loss of generality that the energy £ = 1. We restrict to the case when 
fo G P(R d ) n C(R d ) and we refer to [6] for the relaxation of this condition. 
Since fo G C(M d ), we can define 

N 

with F N (r) := / ]T/ (^)dw, 

J§Nd-l( r ) T" "~ 



so that (ii) holds. 

From the gaussian moment bound M TOo (/q) < oo, we obtain from [22J that there exists 
some constants C > such that 



F N (VN) ~ CN, 
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and for (f(vi, . . . , vg), £ < N: 



Is 



N 



(p{vi, . . . ,v £ ) TT f (vj)dS(V) — >CN / ip(v 1 ,...,ve)dfo(v 1 )...df (ve) 



which proves the chaos. 

We then deduce (i) thanks to |41| Proposition 2.2]. Point (hi) is just a consequence 
of the above asymptotic with the choice ip = mi, t/j = 1. When we assume furthermore 
that /o is compactly supported, say supp/o C [— A, A] d , we deduce supp/^ C {V G 
R Nd , M% 3 (V) < m 3 (A)} and (iv) holds. 

In the non compactly supported case, for any k G N* and for any constant A^ we 
define fo,k '■= fo hv\<A k - It is clear that for any k G N*, there exists f^ k such that f^ k is 
/o, fc -chaotic such that (ii) and (hi) hold and supp/^ C {V G R dN ; M% 3 (V) < m 3 (A k )}. 
For any given sequence (a^) which tends to infinity, we dehne k^ in such a way that 
m 3{Ak N ) = Qn so that kjy — > oo when iV — > oo. We then easily verify that /q := 
satisfies the properties (i)-(iv). □ 

6.2. Chaos propagation without compact support for the Boltzmann equation. 

We may relax the compactly support condition in Theorem 14.11 and Theorem 15.11 thanks 
to Lemma E31 We assume that 



M mo (fo) ■= / e a ^ f (dv) < oo 



for some a G (0, oo) and we dehne /q as in Lemma [6.4i Instead of (|3.4ip in Theorem 13.271 
we have the following bound. For any increasing sequence — > oo, for any ip G (J-\ fl 
J^2 n ^3) , there exists a constant Ci^ (independent of N) such that for any N G N*, 
with N > 21, 



[0,T) 



(6.34) sup (to /( f)-(^(/o)) 87V )^ 

^ + q£U cfe £2 (iV) + e Qjv>r (w^^/o^, 5/o 



By choosing (a^) appropriately, we will deduce from (|6.34p that S^(f^) is again S^ L (fo)- 
chaotic uniformly in time. 

In the case of the Boltzmann model for Hard Spheres or true Maxwell Molecules, we 
can take distg 3 = W\ the usual Monge-Kantorovich-Wasserstein distance in P(R d ). We 
claim that 

(6.35) W Wl (*$f?,5 f0 ) — ► 0. 

iV— >oo 

First, thanks to [41; Proposition 2.2] and the fact that is /o-chaotic from Lemma HT4"1 
we deduce that ftpf^ converges to 5f in the weak sense in P(P(R d )) (that means taking 
duality product with functions of C(P(R d ))). Next, thanks to [Ml Theorem 7.12], ([6T35]) 
boils down to prove that 



(6.36) lim sup [ Wi(pJ o )n%fF(p) = 0, 

R-*oc NeN * JWi(p,fo)>R 

which will be a straightforward consequence of the following boun 

(6.37) sup / [W 1 (^,f )] 2 f^(dV)<oo. 

AfeN* Je n 
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Finally, in order to get (I6.37D . we infer that from |46} Theorem 7.10] 
.N t M2 / w.N t II 2 



< 2\\^f Ml + 2 \\f f Ml 



< 2[Mf(F)] 2 + 2||/ ||^ 1 <2M$ f (V)+2\\f Q f Ml . 



That implies 



[WM,f )] 2 f N (dV) <2\\f \\"+2(f N ,M 2 ), 

E N 1 



which, together with (ii) in Lemma [674], ends the proof of (16. 37ft and then of (|6.35p . 

From the fact that the &t,a functions exhibited in Theorem 14.11 and Theorem 15.11 are 
independent of T and satisfy @t,a{%) = &a(x) — > when x — > for any fixed A G (0, oo) 
we may build (thanks to a diagonal process) a sequence (a ) such that 

e Q iv (vM*£/o^ /o )) — > °- 

iV— >oo 

Coming back to (j6.34f) we obtain that (f^) is S'/ vrL (/o)-chaotic uniformly in time. 
6.3. Chaoticity of the sequence of steady states. 

Theorem 6.5 (Abstract fluctuation estimate in the infinite time). Consider a sequence 
of initial datum which satisfies (Al) (with C^ ms = a N may depend on N) and is fa- 
chaotic with /o G Pq x fl Pg. A . Assume moreover that the assumptions of theorem \3.2T\ hold 
with T = oo. Assume finally that ft j N when t — > oo in the weak sense of measures in 
P(E N ) as well as ft — 7 when t — > oo in the weak sense of measures in P(E). For any 
I G N* and tp G {T\ fl Ti fl ^3) } i/iere exists a constant such that for any N G N*, 
wii/i A > 2f , we have 

(6.38) 



< 



<a 



As a consequence, 7^ is ^-chaotic. 

The proof of that result is trivial: we just have to apply Theorem 13.271 and to pass to 
the limit in the left hand side of the inequality (|3.4ip in order to gat (|6.38j) . Arguing as 
in section 16.21 we deduce 7^ is 7-chaotic (whenever Cq ms = a N grows slowly enough) . 

The application to the Boltzmann equation is the following. Consider a sequence of 
initial data fF such that supp 

jN c §A f rf-i(^/]v) ) an d suc h that $ is / -chaotic with 
/ \ v \ 2 fo = 1) / v i fo = for any i = 1, . . . , d. On the one hand, we know (see [22], E]) that 
f^ converges in the large time asymptotic to 7 , the uniform distribution on the sphere 
§ Nd - 1 (y/N) (that holds in L 2 (j N ) with rate exp(-Ajvt) whenever f* G L 2 (j N )). We also 
know (see [36] and the references therein for the hard sphere case and [421 [7] or section 14.11 
for the (true) Maxwell Molecules case) that ft converges in the large time asymptotic to 
7, the normalized Gaussian. As a consequence, we get ()6.38p . That result may seem to be 
trivial, and it is in some sense, because an explicit computation (which go back at least to 
Mehler in 1866) yields the same result. However, our proof it is not based on an explicit 
computation nor a variationnal/entropy optimization principle. The consequence is that 
it applies to many more situations, in particular in the case of some dissipative Boltzmann 
equation (linked to the Granular media), see 
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6.4. On statistical solutions and the non uniqueness of its steady states. Let us 

consider the TV-particles system associated to the Boltzmann collision process that we do 
not write in dual for as we have done before, it writes 



(6.39) 



N 



i<j 



S d-i 



da. 



We want to describe how the BBGKY (Bogoliubov, Born, Green, Kirkwood and Yvon) 
method introduced to derive Boltzmann's equation from Liouville's equation applies in 
our simpler space homogeneous context. Let us thus also introduce the A;-th marginal: 



ff{vi,...,vi) :-- 



f N (v!, ...,V£,We+i, .-.,Wn) dw e+1 ...dw N 



(N-k) 



Integrating the master equation f|6.39[) leads to 



0{£ 2 /N) 



+ 



— V K z 1 



N 

i<i<j 



N ^ " ij 

i,j>£ 



with 



i(N-e-i) j gN-i 



D 



f N (...,v' i ,...,v' j ,...) - f N dadv e+i ...dv N 



N 



Only the second term does not vanish in the limit TV — > oo, so that assuming that ff -»■ 
tt£, we find that (irg) is a solution to the infinite dimensional system of linear equation 
(the Boltzmann equation for a system of an infinite number of particles or the statistical 
Boltzmann equation) 

(6.40) df e = A e+1 (7r e+1 ) 

with -K£ = TTg(t, v\, ...,vi) > and 



a ■ (vj - v e+ i)) dv e+1 da, 



V GR de ^ A e+1 (ir e+1 )(y) = V / \ n+1 (yj)-* i+1 (y)} b( 

with V- = (v 1 ,...,v' j ,...,v i ,v' i+1 ), v'j = v'(vj,v e+1 ,a), v' e+1 = vl(vj,v i+1 ,a) where vectors 
v' and are defined by (jl.4p . 

Lemma 6.6. There exists a non chaotic stationary solution to the statistical Boltzmann 
equation. In other words, there exists ir G P(P(R d )) such that tt ^ 8 p for some p G P(R d ) 
and Ai + i(tt£ + i) = for any t G N. 

Proof of Lemma 16.61 It is clear that any function on the form 

V g R d( - e+1) i — y tt £+1 (V) = 4>{\V\ 2 ) 

is a stationary solution got equation (|6.4U|) . that is Q+xC^+i) = 0. Now we define, with 
d = 1 for the sake of simplicity, the sequence 



(1 + \y]2^m+i/2 
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with c\ such that n\ is a probability measure and C2 = c\ a 2 with a 2 chosen in the following 
way: 

a 2 , ol 2 f 1 

a«* = — — ^-rr / ^ dv* 



(1 + V 2 + V 2 ) m + 1 (1 + v 2 )' m+l J R Q 1 _^')r?i+l 



1+u 2 

«2 f , 

aw* 



(l+f 2 ) m+1/2 il( 1 + ™*) 2 (1 + «2)m+l/2 ' 

By an iterative process we may chose the constants eg in such a way that tt£ is a solution to 
A-iipi) = (because it is a function of the energy) and satisfies the compatibility condition: 

ne(V)= / n£ + x(y,v*)dv*. 
Jr 

We have exhibit a solution which is not chaotic. □ 
We come back to the abstract setting. We start with the TV-particles system equation, 

d t f N = A N f N , 

that we write in dual form 

Lemma 6.7. Assume that 

(A6) (//) is %fa in P(P^) /or any ^ > N (or equivalently, f N is tight in P(P(E))); 
(A7) Gf +l ip — > Gf +l (p when N — > 00, for any fixed tp G Cb(E e ). 

Then, up to extraction a subsequence, (/ ) converges (in the sense of any l-th marginals) 
to a solution tt = (7^) G P(P(E)) to the infinite Hierachy 

d t 7T = A C °TT in P(P(E)), 

which simply means 

dt{n£, tp) = {TT£ + \,Gf +l Lp) for any I G N*. 
6.5. Uniqueness of statistical solutions and chaos. Assuming (A2i) and 
(A2iii') [0,oo) — > Pg 1 , t >->• S^ L f uniform continuously for any / G Pg 1 , 
which is a weak version of (A2iii), we have that (St) is a co-semigroup. As it has proved 
in step 1 of Lemma 12.131 f° r anv ^ £ Cb(Pg 1 ,M.) we may define T°°<I> by 

(T°°<S>)(f) = $(S? L f), 

and we build in that way a co-semigroup (TJ?°) on (^(Pg^R). The Hille-Yosida theory 
imply that there exists an closed operator G°° with dense domain dom(G°°) in Cb(Pg 1 ,M) 
so that (T t °°) is the semigroup associated to the generator G°° . 

Now, on the one hand, for any ttq G P{Pg 1 ) we may define the semigroup (St) on P(Pg 1 ) 
and the flow (7ft) by setting 7ft = Sf°7TQ and (duality formula) 

V$ G C b (P gi ;R) {SZ°n ,$) = (ir ,T t °°$). 

Remark (see gj) that S t °°7r G (C 6 (P(F))' + ^ P(P(E)). Under the additional as- 
sumption that Pg lt a is compact for any a (that is true in our application cases when 
|| • || g x metrizes the weak measures topology) that relation defined a unique probability 
S^ttq G P(Pg 1 ), and again (S^°) is a co-semigroup on P(Pg 1 ). 

On the other hand, we say that irt G C(R+; P(Pg ± )) is a solution to the equation 

(6.41) d t 7T t = A°° 7T t , 
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if for any $ G dom(G°°) there holds 

^(7r tj *> = fa,G°°*> in V '([0, oo)). 

Theorem 6.8. Assume that (A2) and (A4) hold, as well as that C 1 (Pg 1 ;M) is dense 
in C; ) (Pg 1 ;]R). For any initial datum ttq G P(Pgi) the flow itt is the unique solution in 
C([0, oo); P(Pg 1 )) to (|6,4ip starting from ttq. Moreover, if ttq is fa-chaotic, then ir t is 
St /o- chaotic for any t > 0. 

Proof of Theorem \6.8[ Step 1: Chaos propagation. From Hewitt-Savage's theorem [21], 
for any tt G P(P(E)) there exists a unique sequence (ir) G P(E e ) such that the identities 



JP(E) 



(E) 

^(/)vr(d/) = (vr,^) 



hold for any ip G C^E)® 1 . As a consequence, if 7To is /o-chaotic, 

(4^) = (7f t ,i^) = (7r ,rr<) = (rr<)(/o) 



<(5f L /o) = (Sf L /o,^i) ... (5^ /o, 



which means that 7r| = f? , or equivalently 7Tt = <5/ t , and the statistical solution 7rj is 
/ 4 -chaotic. 

Step 2: Uniqueness. For any £ > and n G N* owe define e := t/n and tk = sk, 
s k = t -t k . Then for any <3? G C^Pg^R) we define $ t ■= T t °°<I>. The very fundamental 
point is that thanks to Lemma 12.131 we have &t 6 C^ ) (Pg 1 ;M,) C dom(C7°°) for any t > 0. 
We write 

(7T t , $) - (7f t , $) = (7T t , $) - (7T , $ t ) 
n-1 

= Z){[< 7r **+l'*'*+l) - >*«*>] + [K+l> $ «*> " *-*)]} 

fc=0 

n-1 

= 7i + r 2 = {7i,k + r 2jfe } . 

fc=0 

On the one hand, we have 

Tx, k = (7r tfc+1 ,$ Sfc+1 -T-$ Sfe+1 ) = -(7r tfc+1 ,^ E ^[Tr$ Sfc+1 ]d S ) 

= [ £ lG°°<S> Sk+1+s ]ds) = - f k (7r t _ [s+lj£] ,G°°$ s )d S , 



where [s,e] = [s/e] e. Passing to the limit n — > oo, we get 

*t ft 



71 = - / (7r t _ [s+1)£] ,G°°$ s )d S — > - / <7ri_ fl ,G°°$ fl )ds. 
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On the other hand, we have 

f £ d 

T 2 ,k = ^ —{TT tk+T ,$ Sk )dT 

= [ E (Tr tk+T ,G°°$ Sk )dT 
Jo 

= / fc+1 (7r T ,G°°^_ M )dr. 

Passing to the limit n — > oo, we get 

T 2 = f t (7r r ,G °$ t _ [r)e] )dr — > f\^G°°^ T )dT. 
Jo n ^°° Jo 

As a conclusion, for any $ £ C 1 (Pg 1 ;IR), we have proved 

(7T tj $) = (7f t ,*). 

From a density argument we conclude that lit = ^t- D 
Gathering Lemma 16.71 and Theorem 16.81 we obtain a propagation to the chaos result. 

Corollary 6.9 (Abstract chaos propagation). Assume (A2), (A4) ; (A6), (A7) and the 

following compatibility between G°° and the sequence G^°: 

VfGN*, V<peJf e (G N ,R V ) — > (G°°,R V ). 
Assume furthermore that f^ is fo-chaotic. Then f^ is Sj^ L fo-chaotic. 
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